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ABSTRACT 


The  compressive  flow  of  viscoelastic  materials 
between  two  parallel  flat  disks  under  a  constant  load  has 
been  investigated  analytically,  numerically,  and  experi¬ 
mentally.  This  process  simulates  a  number  of  compression 
molding  and  lubrication  experiments;  the  purpose  of  our 
study  was  to  assess  the  effects  of  fluid  viscoelasticity  and 
of  temperature  .gradients  in  these  applications. 

A  dimensionless  group  3  {=  )  has  been  found 

/mih  R 

very  useful  for  determining  the  flow  regimes  when  there 
exists  a  substantial  transverse  viscosity  gradient  in  the 
fluid  charge,  such  as  in  the  nonisothermal  compression 
molding  processes. 

Compressive  flow  of  linear  viscoelastic  materials 
has  been  analyzed  analytically.  It  shows  that  the  squeezing 
motion  becomes  oscillatory  when  the  ratio  of  the  Deborah 
number  to  the  Reynolds  number  is  larger  than  a  critical 
value,  and  that  tie  linear  viscoelastic  materials  are 
squeezed  faster  than  the  cor.’*ssponding  Newtonian  cases. 

Compressive  flow  of  various  non-linear  model  fluids 
has  also  been  analyzed  numerically.  The  Maxwell  fluid 
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XV 


behaves  much  like  linear  viscoelastic  materials,  except 
under  extraordinarily  high  loading  conditions.  But,  the 
Johnson-  Segalman  model  and  the  Marrucci  structural  model 
show  that  slower  squeezing  may  arise  after  the  initial  rapid 
transient  under  moderate  loading  conditions.  This  slower 
squeezing  must  be  due  to  the  special  features  of  these 
models,  which  the  Maxwell  model  does  not  exhibit,  such  as 
stress  overshoot  in  the  transient  flows. 

Experimentally  two  different  observations  have  been 
made.  A  silicone  polymer  shows  the  oscillatory  and  the 
faster  squeezing,  which  is  predictable  by  the  Maxwell  type 
of  model  fluid.  Two  other  polymer  solutions  show  an 
inflection  point,  which  probably  reflects  a  very  weak 
oscillation,  and  a  slower  squeezing  than  for  the 
corresponding  inelastic  cases.  The  slower  squeezing  of 
these  solutions  seems  to  be  due  to  the  transient  behavior  of 
these  materials  such  as  stress  overshoot.  The  use  of  those 
models,  which  can  predict  the  transient  behavior  more 
precisely,  is  recommended  to  describe  the  transient 
responses  of  these  materials. 
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CHAPTER  1 


INTRODUCTION 

1 . 1  Relevance  of  the  problem 

The  problem  to  be  considered  in  this  work  is  that  of 
the  compressive  flow  of  viscoelastic  fluids  between  two 
horizontal  circular  flat  disks,  shown  schematically  in 
Fig.  1.1.  .The  test  fluid  is  contained  between  two  disks 
which  are  at  rest  for  times  t<07  at  t=0  the  upper  disk  is 
released  and  falls  under  the  normal  load  F.  The  spacing 
between  the-  disks  is  measured  as  a  function  of  time. 

This  compressive  flow  between  disks  is  of  interest 
for  many  reasons . 

(1)  It  is  encountered  in  the  popular  plastometer 
method,  which  has  been  used  to  determine  the  material 
properties  of  highly-viscous  materials. 

(2)  It  is  also  encountered  in  certain  polymer 


*The  terra  “compressive  flow"  stands  for  the  flow  in  the 
opposite  sense  of  the  extensional  flow,  not  for  the  flow  of 
compressible  fluids.  It  is  equivalent  to  "squeezing  flow" 
or  "squeeze  film  flow"  which  has  been  used  more  often  in  the 
literature . 
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F 
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Fig.  1.1  Schematic  diagram  of 
the  compressive  flow. 
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processing  operations  such  as  compression  molding,  injection 
molding,  stamping,  etc.  The  polymeric  charges  in  these 
processes  are  frequently  filled  with  fibers,  whose 
orientation  will  determine  the  material  properties  of  final 
products,  and  the  orientation  of  the  fibers  is  believed  to 
be  determined  by  the  flow  behav.ior  of  the  polymeric  medium. 

(3)  It  also  arises  in  rubrication  systems,  and  there 
has  been  a  controversy  as  to  whether  or  not  viscoelastic 
lubricants  will  perform  better  than  Newtonian  lubricants, 
since  modern  motor  oils  often  contain  polymeric  additives 
which  render  them  viscoelastic.. 

(4)  Most  of  all,  this  flow  is  of  particular  interest 
to  rheologists  since  both  shearing  and  extensional 
deformations  are  present  under  transient  conditions,  the 
flow  being  dominated  by  shear  near  the  wall  and  by  extension 
in  the  middle  of  the  gap.  Therefore,  the  compressive  flow 
is  a  good  candidate  to  be  used  in  evaluating  proposed 
constitutive  equations,  especially  transient  responses,  and 
further  improving  them. 

It  will  be  assumed  throughout  that  the  fluid  is 
incompressible . 

1 . 2  Review  of  previous  work 


In  this  section  we  will  review  the  principal 
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theoretical  studies  and  experimental  work  which  have  been 
published . 

Stefan(1874)  appears  to  be  the  first  to  have  dealt 
with  this  problem.  He  analyzed  the  flow  for  Newtonian 
fluids  and  derived  the  equation 

j_ 

2 

(1.1) 


H(t)  = 


Ho^ 


I6F  t 
3717 


which  is  known  as  the  Stefan  equation  after  him.  A  complete 
list  of  symbols  and  their  meanings  are  found  in  the 
Nomenclature.  A  more  systematic  derivation  of  equation 
(1.1),  based  upon  the  creeping  flow  approximation  and  the 
parallel  squeezing  assumption,  is  given  by  Denn(19a0).  Here 
the  term  "parallel  squeezing"  implies  that  material  planes 
which  are  initially  horizontal  remain  so  during  the 
subsequent  deformation.  This  assumption  is  not  always 
valid,  even  though  it  is  very  useful  in  many  cases.  In 
Section  2.2,  we  will  consider  an  example  in  which  the 
parallel  squeezing  assumption  fails.  This  has  been  also 
mentioned  by  Brindley  et.  al.(1976).  The  Stefan  equation 
has  been  tested  and  verified  experimentally  for  Newtonian 
fluids  by  many  authors,  including  Parlato( 1969 ) , 
Leider ( 1974) ,  Brindley  et.  al.(1976),  and  Griram( 1977 ) . 

In  1931,  the  squeezing  of  power-law  fluids  was 
analyzed  by  Scott(1931),  who  developed  an  equation  similar 


to  equation  (1.1)  given  by 
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H(t)  = 


n+i 

XI  "  4.  2(n-H)  F(nt3)  f 
Ho/  .^n  +  i  71  K  ^ 


n 

n+ 1 


(1.2) 


where  K  and  n  are  the  consistency  factor  and  the  power-law 
index  of  the  given  power-law  fluid,  respectively.  The 
derivation  of  equation  (1.2)  will  not  be  repeated  here  since 
complete  developments  of  the  Scott  equation  are  given 
elsewhere  (Scott , 1931 ;  Leider  and  Bird, 1974;  Grimm, 1977). 
Experimental  work  on  power-law  fluids  has  been  carried  out 
by  Parlato(i969) ,  Leider(1974) ,  Brindley  et .  al.(1976),  and 
Grimra(1977).  In  general,  experimental  data  agree  well  with 
the  Scott  equation. 

Various  researchers  have  considered  the  case  of 
viscoelastic  fluids  in  the  compressive  flow.  Ail  of  them 
except  Metzner( 1968a)  predict  that  viscoelastic  fluids 
squeeze  out  faster  than  the  corresponding  inelastic  fluids, 
which  is  the  opposite  of  many  available  experimental 
results . 

Tanner (1965)  analyzed  the  flow  for  contravariant 
convected  Maxwell  fluids  with  a  power-law  viscosity  and  a 
constant  relaxation  time.  He  argued  that  the  normal  stress 
effects  are  small  compared  to  the  shear  stress  effects  and 
predicted  faster  squeezing  of  viscoelastic  fluids  than  the 
corresponding  power-law  fluids. 
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Metzner( L963a)  appears  to  be  the  first  to  have 
recognized  the  possible  importance  of  the  extensional  flow 
as  well  as  the  shearing  flow  in  tljis  problem.  He  mentioned 
that  extraordinarily  high  stresses  are  predicted  to  be 
reguired  for  rapid  extensional  defoirmation  of  viscoelastic 
fluids,  and  he  predicted  slower  squeezing  of  viscoelastic 
fluids  ( contravariant  convacted  Maxwell  fluids)  based  upon 
the  "extensional  primary  field"  approximation 
(Metzner , 1971 ) .  Williams  and  Tanner(1970)  also  considered  a 
combination  of  shear  and  extensional  effects  but  concluded 
that  extensional  effects  are  small  compared  to  the  shearing 
effects . 

Kramer’s  analysis (1974 )  is  unique  in  that  the 
particle  path  equations  are  numerically  solved  using  a 
convectsd  coordinate  system  without  neglecting  any  of  the 
normal  stresses.  He  used  the  integral  constitutive  equation 
of  Lodge's  rubberlike  liquid(see  Lodge,  1964)  with  a  single 
exponential  memory  function,  which  is  identical  to  the 
contravariant  convectsd  Maxwell  fluid.  By  assuming  parallel 
squeezing  and  negligible  inertia,  he  predicted  an  initial 
instantaneous  drop  and  more  rapid  squeezing  of  viscoelastic 
fluids.  The  initial  drop  implies  that  the  material  has  no 
resistance  at  the  instant  t=0  except  its  own  inertia;  no 
resistance  of  the  material  implies  zero  apparent  viscosity, 
hence  infinite  Reynolds  number.  Thus,  the  inertia  should  be 
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taken  into  account  in  the  compressive  flow  of  viscoelastic 
materials  under  a  constant  load,  even  though  the  creeping 
flow  approximation  (negligible  inertia)  is  very  useful  in 
the  inelastic  cases. 

Leider  and  Bird(1974)  have  suggested  that  the  use  of 
a  rheological  model  which  can  describe  stress  overshoot  in 
simple  shear  flow  is  imperative  to  explain  the  slower 
squeezing  of  viscoelastic  fluids.  The  models  used  in  the 
three  analyses  cited  above  do  predict  the  presence  of  a 
first  normal  stress  difference,  but  do  not  predict  the 
stress  overshoot  phenomena  or  a  second  normal  stress 
difference  in  shear  flows .  Leider  and  Bird  proposed  an 
empirical  equation  with  an  overshoot  correction  factor  to  be 
used  in  squeeze  film  problems. 

Brindley  at.  al.(1976)  analyzed  the  flow  for  the 
second  order  fluid  and  again  predicted  faster  squeezing  of 
viscoelastic  fluids.  The  second  order  fluid  shows  the  first 
and  second  normal  stress  differences  in  shear  flows,  but  no 
stress  overshoot.  They  also  presented  some  very  interesting 
experimental  results  which  show  solid-like  bouncing  behavior 
under  some  severe  loading  conditions,  but  they  did  not 
provide  a  theoretical  explanation  of  this  bouncing  behavior. 

Experimental  studies  on  viscoelastic  materials 
include  those  of  Parlato( 1969 ) ,  Leider ( 1974) ,  Brindley 
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et.  al.(1976),  and  Grimra(i97a ) .  In  general,  with  slow 
squeezing  (low  Deborah  number),  the  squeezing  behavior  of 
viscoelastic  fluids  is  close  to  the  corresponding  inelastic 
case.  With  fast  squeezing  (high  Deborah  number)  most  of 
experimental  results  show  that  viscoelastic  materials 
squeeze  out  much  slower  than  the  corresponding  inelastic 
fluids,  and  under  very  severe  loading  conditions  some 
materials  even  bounce  back  after  some  amount  of  squeezing. 

Some  authors  (Leider , 1974;  Leider  and  Bird, 1974) 
have  used  the  "half  time",  which  is  the  time  required  for 
the  disks  to  move  from  a  separation  Ho  to  Ho/2,  to  represent 
their  experimental  results.  The  "half  time"  is  very  useful 
in  representing  the  data  of  purely  viscous  materials,  but  it 
is  not  recommended  in  the  viscoelastic  cases  since  it  may 
conceal  interesting  elastic  effects.  This  has  been  pointed 
out  by  Binding  et.  al. (1976a). 

Tichy  and  Winer (1978)  studied  constant  speed 
squeezing  flow  instead  of  constant  load  squeezing,  using 
Lodge's  rubberlike  liquid  model  with  a  single  relaxation 
time.  They  predicted  that  the  load-bearing  capacity  of  a 
viscoelastic  fluid  may  be  increased  due  to  normal  stress 
effects  or  decreased  due  to  a  delayed  response  of  shear 
stress  to  a  change  in  shear  rate. 

Shirodkar ( 1981 ) ,  and  Shirodkar  and  Middleman ( 1982 ) 


also  considered  constant  speed  squeezing  flow.  They  used  a 
constitutive  equation  due  to  Wagner ( 1976 ) ,  which  is  an 
empirical  modification  to  the  Lodge's  rubberlike  liquid 
model.  Even  though  Wagner's  model  is  capable  of  exhibiting 
non-Wewtonian  viscosity,  normal  stress  in  simple  shear  flow, 
and  stress  overshoot  in  transient  simple  shear  flow,  the 
principal  drawback  in  this  model  is  the  absence  of  a  general 
form  of  the  damping  function.  They  predicted  that  fluid 
elasticity  can  increase  the  force  resisting  the  approach  of 
the  boundaries  of  a  squeeze  film  at  high  shear  rate.  They 
also  suggested  that  generalizations  regarding  the  role  of 
viscoelasticity  may  be  impossible  and  the  effect  of 
polymeric  additives  on  load-bearing  capacity  appears  to 
depend  upon  whether  the  motion  is  under  constant  load  or 
constant  speed,  or  some  combination  thereof. 

The  properties  of  polymeric  materials  responsible 
for  the  slower  squeezing  found  experimentally  are  still  in 
question  and  need  to  be  studied  more  systematically. 

There  has  been  a  somewhat  different  question  in  the 
compression  molding  process  (usually  nonisothermal)  which 
concerns  the  flow  in  the  mold.  It  is  important  to  identify 
the  portion  of  the  material  that  flows  preferentially  to 
fill  the  remainder  of  the  mold  cavity:  the  fluid  near  the 
wall  or  that  in  the  central  portion  of  the  mold.  The  answer 
will  depend  upon  factors  such  as  geometry,  material 
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properties,  and  temperature  gradient  in  the  mold.  Both 
cases  have  been  observed  experimentally.  There  has  been  no 
quantitative  analysis  of  this  flow. 

1.3  Objectives  and  approaches  of  present  work 

Several  basic  problems  in  compressive  flow  remain 
unsolved.  Resolving  those  problems  will  be  the  main 
objectives  of  this  work,  as  follows: 

(1)  We  require  a  better  qualitative 
understanding  of  the  behavior  in  compressive 
flow  of  viscoelastic  materials  between  disks 
under  isothermal  conditions. 

(2)  Furthermore,  we  wish  to  analyze  this  flow 
quantitatively  to  explain  the  experimental 
observations  of  slower  squeezing  and  bouncing 
behavior  of  viscoelastic  materials. 


(3)  Finally, 

in 

the 

compression 

molding 

processes , 

we 

wish 

to 

understand 

the 

preferential 

flow 

of 

some 

portion 

of 

the 

material  in  the  mold  cavity. 

In  chapter  2  we  will  consider  the  compressive  flow 
of  purely  viscous  fluids,  including  Newtonian  and  power— law 
fluids.  The  flow  in  the  mold  cavity  will  be  also  examined 
quantitatively  in  this  chapter. 
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Compressive  flow  between  disks  is  not  simple;  it  is 
a  combination  of  shearing  and  extensional  flow  under 
transient  conditions,  and  the  viscoelastic  nature  of  the 
materials  of  primary  interest  makes  the  problem  even  more 
complex.  Before  dealing  with  this  complex  problem  we 
consider  a  rather  simple,  but  closely  related  problem  in 
Chapter  3;  this  is  lubricated  compressive  flow.  Since 
lubricated  compressive  flow  between  disks  is  essentially 
biaxial  extensional  flow,  various  kinds  of  non-linear 
viscoelastic  constitutive  equations  can  be  treated  without 
great  difficulty.  The  unlubricated  problem  will  be 
considered  in  Chapter  4.  Numerical  computation  is  required 
in  this  case.  The  finite  element  numerical  technique,  which 
will  be  tested  in  the  simple  lubricated  problem  in  Chapter 
3,  will  be  applied  to  solve  the  continuity,  momentum,  and 
constitutive  equations  simultaneously.  In  Chapter  5, 
experimental  results  will  be  compared  to  the  theoretical 
predictions . 


CHAPTER  2 


COMPRESSIVE  FLOW  OF 
PURELY  VISCOUS  FLUIDS 

The  compressive  flow  of  purely  viscous  fluids  will 
be  treated  in  this  chapter,  focusing  on  the  following: 

(1)  We  will  investigate  the  validity  of  the 
Stefan  equation  (for  a  Newtonian  fluid)  and 
Scott  equation  (for  a  power-law  fluid)  by 
finite  element  calculations,  in  Section  2.1  and 
2.3,  respectively. 

(2)  We  will  analyze  the  flow  in  the  mold  cavity 
for  the  compression  molding  process  in  Sections 
2.2  and  2.4. 

(3)  Throughout  this  chapter  we  will  test  the 
finite  element  and  other  numerical  schemes  used 
to  solve  compressive  flow  problems. 

2. 1  Newtonian  fluid 

The  compressive  flow  of  an  isothermal  inelastic 
fluid,  neglecting  fluid  inertia,  is  governed  by  the 
following  equations  and  boundary  conditions: 
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Continuity ; 


Navisr-Stokes  (variable  viscosity,  7  =  |(r..)): 


0  = 

ar 

(2.2) 

0  = 

9Z 

(2.3) 

Boundary 

conditions  (see 

Fig.  2.1  for 

description 

of  the 

coordinate 

system  used) : 

Vr  = 

0  , 

o 

II 

at  r  =  0 

(2.4) 

Vz  = 

0  , 

=  0 
dz 

at  z  =  0 

(2.5) 

II 

0  , 

II 

1 

< 

at  z  =  H 

(2.6) 

All  the  stresses  vanish  at  the  free  (2.7) 

surface 

The  second  condition  in  (2.4)  comes  from  the 
requirement  of  zero  shear  stress  along  the  axis  of  symmetry. 
In  the  case  of  the  parallel  squeezing  assumption (Section 
2.2.1)  this  condition  is  automatically  satisfied.  This 
problem  has  been  solved,  based  upon  the  parallel  squeezing 
assumption,  to  give  the  Stefan  equation  (1.1). 


constant  force 


Fig.  2.1  The  domain  used  in  the 
numerical  calculations . 
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We  wish  to  solve  the  problem  niamerically  with  no 
assumptions  in  order  to  compare  the  solution  to  the  Stefan 
equation.  A  finite  element  technique  will  be  used:  the 
finite  element  formulation  for  generalized  Newtonian  fluids 
is  given  in  appendix  D.l. 

Because  of  flow  symmetry,  the  domain  of  interest  in 
numerical  calculations  will  be  confined  to  a  quarter  of 
entire  domain  (shaded  portion  in  Pig.  2.1).  Neglecting 
fluid  inertia,  the  problem  becomes  a  quasi-steady  state 
problem  even  though  the  flow  itself  is  transient. 
Therefore,  at  each  time  instant  we  solve  the  steady  state 
equations  (2. 1-2. 3)  under  the  boundary  conditions  (2.4-2. 7) 
on  the  given  domain.  Then,  after  a  given  time 
increment(At  ),  we  move  the  boundary  of  the  domain  by 
V  *  At  »  solve  the  problem  on  the  new  domain,  and  proceed  to 
the  next  time  step.  The  numerical  algorithm  is  given  in 
Fig.  2.2.  Between  time  steps,  a  predictor-corrector  method 
is  used;  since  the  problem  is  entirely  linear,  the  top 
plate  velocity(V)  is  adjusted  linearly  to  produce  the  given 
constant  force. 

Numerical  calculations  have  been  carried  out  for 
three  different  values  of  R/H© (=5, 15, 50) .  The  results  are 
compared  to  the  Stefan  equation  in  Fig.  2.3,  in  which 
dimensionless  film  thickness (H/H© )  is  plotted  against 
dimensionless  time,  t(= -^(-^jy)  .  At  large  values  of 


Read  data  (  X(I),  FORCED,  ,  etc.) 

I 

Time  zero  ?  -S2 - 

jyes 

Finite  element  solution  (V(I),  p(I)) 

I  ,  ~ 

Calculation  of  force  (FORCE) 


Linear  adjustment  of  solution 
y  ( I )  =\J  (I )  *  (  FORCEO/ FORCE ) 
p(I)=p{I)*(FORCEO/FORCE) 


Print  solution 


Next  time  step  ? 


■STOP 


Sold(^)=2(l) 


Read  previous 
solution  (V(I),p(l)) 


Enter  ^t 


Move  boundary  nodes  (predictor  part) 
Xnew  *'Xoid+:^oid*^t 


Finite  element  solution  (Vj^^^d), 

1 

Calculation  of  force  (FORCE) 


Linear  adjustment  of  solution 


Move  boundary  nodes  (corrector  part) 
Snew=^2old+i<Vold 


Fig.  2.2  Numerical  algorithm  to  solve  compres¬ 
sive  flow  of  Newtonian  fluids. 


Pig.  2.3  Dimensxonless  plate  spacing  versus 
dimensionless  time  for'  isothermal  squeezing 
of  Newtonian  fluids. 


R/Hq ,  numerical  results  agree  well  with  the  Stefan  equation, 
but  as  R/Hq  decreases  the  discrepancy  between  the  numerical 
solution  and  the  Stefan  equation  increases.  This 
discrepancy  is  due  to  the  velocity  rearrangement  caused  by 
the  stress  singularity  at  the  edge  of  the  disk  and  becomes 
small  as  R/H©  increases,  which  is  also  seen  in  the  pressure 
profiles  (see  Pig.  6  and  7  in  Appendix  A).  In  other  words, 
the  parallel  squeezing  assumption  is  a  good  one  as  long  as 
R/H  is  large  enough. 

In  the  next  section  we  will  see  a  typical  example  in 
which  the  parallel  squeezing  assumption  no  longer  holds  and 
the  boundary  condition  at  the  edge  of  the  disk  plays  an 
important  role . 

2.2  Compr  e  s  s i v  e  flow  between  parallel  disks  with  a 
transverse  viscosity  gradient 

Let  us  consider  the  flow  in  the  mold  cavity  in  the 
compression  molding  process,  which  is  depicted  schematically 
in  Fig.  2.4.  There  have  been  two  different  observations  on 
the  flow  patterns  induced  as  the  mold  is  closed;  Marker  and 
Ford (1977)  observed  the  preferential  flow  of  hot  material 
near  the  mold  faces,  while  Denton ( 193.’. )  and  Denn,  Tadraor, 
and  Edelist ( 1981 )  reported  no  preferential  flow  of  the  fluid 
near  the  wall  and  found  the  maximum  velocity  at  the  center 
plane.  These  flow  patterns  are  shown  schematically  in 
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Fig.  2.4  Schematic  diagram  of  the 
compression  molding  process. 
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Fig.  2.5.  Qualitatively  it  is  obvious  that  both  cases  are 
possible,  depending  upon  the  material  properties,  geometry, 
and  processing  conditions. 

A  detailed  quantitative  analysis  of  this  problem  is 
given  in  Appendix  A,  and  only  its  essence  will  be  discussed 
in  this  section.  Here  we  consider  the  flow  between 
approaching  disks  when  there  is  a  substantial  viscosity 
gradient  in  the  fluid  charge;  this  viscosity  difference  may 
result  from  the  temperature  difference  caused  by  the  hot 
plates,  or  it  may  represent  an  approximation  to  the 
properties  of  a  viscoelastic  charge  in  which  the  resistance 
to  a  biaxial  deformation  will  be  much  greater  than  to 
shearing . 

2.2.1  Parallel  squeezing  assumption 

As  already  mentioned.  The  parallel  squeezing 
assumption  has  been  a  conventional  approach  in  the  squeeze 
film  problem.  It  can  be  shown  that  it  is  a  direct 
consequence  of  the  parallel  squeezing  assumption  that  low 
viscosity  fluid  near  the  disks  cannot  flow  out 
preferentially,  regardless  of  the  viscosity  difference 
between  the  center  plane  and  the  wall  (see  Appendex  A). 
Therefore,  we  can  see  that  the  parallel  squeezing  assumption 
will  not  be  valid  when  the  viscosity  difference  and  geometry 
are  such  as  to  generate  a  preferential  flow  of  the  low 


Fig.  2.5  Flow  patterns  observed 
experimentally. 

(a)  preferential  flow  of  centra.,  fluid. 

(b)  preferential  flow  of  the  fluid 
near  the  wall. 
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viscosity  fluid  near  the  wall. 

2.2.2  Squeezing  force  and  boundary  condition 


The  force  required  to  bring  the  plates  together  is 
obtained  by  integrating  the  axial  total  stress  over  the 
surfaces  of  the  disks.  The  calculation  of  this  force 
requires  a  knowledge  of  the  boundary  condition  at  the  outer 
edge  of  the  disks. 


The  usual  boundary  condition  employed  at  the  disk 
edge  is  shown  schematically  in  Fig.  2.6(a);  the  normal 
stress  at  the  free  surface,  (3^2=-p+'<»  ,  is  balanced  against 
atmospheric  pressure.  This  gives 


(2.8) 


For  the  case  of  constant  viscosity,  we  obtain 


(2.10) 


At  constant  force,  this  leads  to  the  Stefan  equation . 


This  boundary  condition  gives  a  paradoxical 
prediction  for  the  case  in  which  there  is  a  large  viscosity 
difference  between  the  central  plane  and  the  walls.  From 


(2.8)  and  (2.9)  the  force  required  is  of  the  order  of  the 
smallest  viscosity  that  exists  over  the  finite  portion  of 
the  gap,  even  if  the  central  portion  of  the  gap  contains 
fluid  of  arbitrarily  large  viscosity.  Obviously  this  can 
not  be  true . 

A  better  boundary  condition  at  the  outer  edge  is  the 
vanishing  of  all  stresses  on  the  free  surface  of  the 
extruded  sheet,  shown  schematically  in  Fig.  2.6(b). 
Assuming  that  the  velocity  rearrangement  is  restricted  to  a 
small  neighborhood  near  the  edge,  we  can  approximate  this 
condition  by  requiring  that  the  net  radial  force  component 
be  zero  at  r=R;  that  is, 

/H 

dz  =  0  at  r  =  R  (2.11) 

J  0 

With  the  parallel  squeezing  assumption,  this  gives 


The  second  term  in  (2.12)  corresponds  to  the  stress  from  the 
biaxial  extension  of  the  central  viscous  fluid. 


When  viscosity  is  constant. 


7LC, 


(2.13) 


As  long  as  ^<<1,  (2. 10) and  (2.13)  are  identical.  Thus,  the 


Stefan  equation  is  unchanged  by  the  alternate  boundary 
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condition . 


The  situation  is  quite  different,  however,  when 
there  is  a  large  viscosity  variation  across  the  gap.  The 
first  term  in  (2.12)  is  of  order  and  the  second  is 

%  W 

of  order  - .  Thus,  the  first  or  second  terra  dominates 

depending  upon  whether  the  group 


'h  max 
7  min 


(2.14) 


is  small  or  large,  respectively,  compared  to  unity. 


2.2.3  The  dimensionless  group  S 

When  S  is  small  coitipared  to  unity,  the  stress  from 
the  biaxial  extension  of  the  center  fluid  is  negligible  and 
the  maximum  velocity  occurs  at  the  center  plane.  Thus  the 
parallel  squeezing  assumption  and  the  conventional  boundary 
condition  are  both  valid. 


As  S  increases,  the  new  boundary  condition  is 
necessary  to  compute  the  correct  order  of  magnitude  for  the 
force,  although  the  parallel  squeezing  assumption  is  still 
valid . 


When  S  is  large  compared  to  unity,  the  parallel 
squeezing  assumption  breaks  down  since  the  maximum  velocity 
occurs  in  the  low  viscosity  fluid  near  the  disks.  This  has 
been  shown  numerically  for  the  case  of  two  fluids;  i.e.,  in 
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which  the  viscosity  is  constant  in  each  of  two  regions,  but 
changes  discontinuousLy  at  an  interface  that  is  initially 
planar  (see  appendix  A  for  the.  details).  Thus,  a  new 

analysis  is  necessary  in  this  case.  In  the  two  fluids  case, 

new  assumptions  other  than  the  parallel  squeezing  assumption 
have  been  made  to  derive  an  analytical  solution  which  is 
found  to  be  in  good  agreement  with  the  finite  element 
numerical  solution. 

Thus,  we  now  have  a  quantitative  criterion  for  the 
two  different  flow  regimes  observed  experimentally;  that 
is, 

*  S  <<  1  flow  regime  of  Fig.  2.5(a) 

*  S  >>  1  ;  flow  regime  of  Fig.  2.5(b) 

2 . 3  Power-law  fluid 

In  this  section,  the  Scott  equation  (for  a  power-law 
fluid)  will  be  tested  by  numerical  computation,  just  as  the 
Stefan  equation  has  been  tested  in  section  2.1. 

Governing  equations  and  boundary  conditions  are  the 
same  as  in  Newtonian  case,  except  that  the  viscosity  is 
given  by 

n-i 

=  K|4-Id|  ^  (2.15) 


where 
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Id  =  -  x[(t>-d)'*  -  (2.16) 

and 

/  <^rr  0  dj-z  > 

d  =  0  dee  0  (2.17) 

\  d  rz  0  dzz  / 

The  fluid  is  assumed  to  be  incompressible,  so  that 

Ij  =  <1„‘  (2-13) 

In  solving  the  problem  numerically,  an  iterative 
scheme  is  used  to  evaluate  the  viscosity  function.  That  is, 
the  viscosity  is  initially  based  on  the  previous  solution 
for  an  earlier  time  step.  One  then  solves  for  the  new 
velocity  and  pressure  profiles,  calculates  the  viscosity 
from  this  new  velocity  field,  and  repeats  this  procedure 
until  the  solution  (velocity  and  pressure  field)  converges 
within  a  given  error  allowance.  The  numerical  algorithm  is 
given  in  Fig.  2.7. 

Numerical  computations  have  been  carried  out  for  two 
different  values  of  R/Ho(5,15)  and  two  different  values  of 
the  power-law  index ( 0 . 7, 0. 5 ) .  The  results  are  given  in 
Fig.  2.8,  in  which  the  Newtonian  case  (n=1.0)  is  also 
plotted  for  the  comparison.  In  the  figure, 
dimensionless  time  t  is  defined  by 


the 
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Read  data  (  X(I),  FORCED,  ,  etc.) 

I  ... 


Time  zero  ?  _ 


Finite  element  solution  ,  pd))" 

^Calculation  of  force  (FORCE) 


Adjustment  of  solution 
V (I ) =V ( I ) * ( FORCEO/FORCE ) 
p(l)=p(l)*(FORCEO/FORCE) 


I 

Convergence  test  ? 

|y®s 

Print  solution 


Next  time  step  ?- - ^STOP 

lyes 

y  y 

(I)^X(I)  •  Read  previous 

Void  ^  ^  ^ ^  ^  solution  (V(l),p(I)) 

i 


Enter  At 


Move  boundary  nodes  (predictor  part) 
~new  ~~old  "^^old  * 


Finite  element  solution  (V„^(I), 
Calculation  of  force  (FORCE; 


Adjustment  of  solution 

no 

Convergence  test  ?  - 

jyes 

-Move  boundary  nodes  (corrector  part) 


Fig.  2.7  Numerical  algorithm  to  solve  compres¬ 
sive  flow  of  power-law  fluids. 


A  ® 


R/Ho=5 


2n 

with  these  definitions,  the  Scott  equation  as  well  as  the 
Stefan  equation  become 

i  =  (  I  t  t)‘^  (2-21) 

no 

Again  we  can  see  that  the  edge  effect  is  more 
significant  at  smaller  values  of  the  R/H©  ratio.  The  effect 
increases  as  the  power-law  index  decreases,  or  as  shear 
thinning  behavior  increases.  This  can  be  explained  by 
considering  the  viscosity  difference  between  the  fluids 
inside  and  outside  the  edge  of  the  disk.  As  the  power-law 
index  decreases,  the  fluid  outside  the  edge  has  a  higher 
viscosity  than  the  fluid  inside  the  edge,  which  implies  that 
the  fluid  inside  experiences  more  resistance  from  the  fluid 
outside.  Therefore,  the  squeezing  speed  becomes  slower  as 
the  power-law  index  decreases.  This  effect  decreases  as 
R/Hc  increases,  and  at  large  values  of  R/H©  the  Scott 
equation  is  still  adequate. 

2.4  Partially-filled  compressive  flow  of  Newtonian  fluids . 

In  this  section  we  consider  the  flow  with  a  moving 
front  in  the  mold  cavity,  depicted  in  Fig.  2.9.  We  will 
focus  on  the  moving  free  surface  and  the  flow  near  this 


(2.19) 

(2.20) 
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at  time  =  0. 


1 


Fig.  2.9  Schematic  diagram  of  the 
flow  of  a  Newtonian  fluid  in  the 
mold  cavity  under  the  isothermal 
condition . 
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moving  front.  Numerical  calculations  are  therefore 
inevitable,  and  the  proper  numerical  scheme  has  to  be 
developed  to  determine  the  new  contact  point  between  the 
fluid  and  the  surface,  as  well  as  the  shape  of  the  free 
surface,  for  each  time  step.  It  will  be  assumed  that  the 
fluid  is  Newtonian  and  isothermal. 

2.4.1  Numerical  scheme 

At  any  instant  of  time,  the  velocity  and  pressure 
distributions  can  be  obtained  by  finite  element  calculation, 
as  long  as  the  domain  is  given  beforehand  At  time  t=0,  we 
start  with  a  domain  of  rectangle  shape  with  a  flat  free 
surface,  as  shown  in  Fig.  2.9,  and  solve  for  the  velocity 
and  the  pressure. 

After  a  small  time  increment ( it ) ,  the  boundary  of 
the  domain  is  changed  by 

^new  ~  (2.22) 

The  new  contact  point  is  determined  by  quadratic 
interpolation  of  three  adjacent  nodes,  which  is  shown  in 
Pig.  2.10.  That  is,  after  At  the  nodal  point  1  moves  to  1^, 
2  to  2',  and  3  to  3^,  etc.,  and  a  point  p  between  1  and  2  will 
move  and  just  contact  the  wall;  point  p^ is  the  new  contact 
point  between  the  fluid  surface  and  the  solid  wall.  The 
solution  and  the  boundary  shape  are  then  improved  through 
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Fig.  2.10  Movement  of  the  frontal 
nodes  and  contact  point. 

... - front 

-  new  front  after  At 
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the  predictor-corrector  method.  These  procedures  are 
repeated  at  each  time  step.  This  numerical  algorithm  is 
given  in  Fig.  2.11. 

2.4.2  HJumerical  calculation  results 

Two  sets  of  calculations  have  been  carried  out,  one 
at  a  small  value  of  Ro/Ho(=5),  the  other  at  a  large  value  of 
Ro/Ho (=15) . 

In  the  case  of  Ro/Ho=5,  the  movement  of  the  boundary 
and  the  frontal  free  surface  are  shown  in  Fig.  2.12  at 
various  stages  of  squeezing.  The  radial  velocity  profiles 
of  the  present  case  are  compared  in  figures  2.13  through 
2.16  to  the  corresponding  fully-filled  case,  which  has  an 
extra  amount  of  fluid  outside  the  edge  of  the  disks;  the 
radius  for  the  fully-filled  calculation  is  based  on  the 
contact  point,  p.  At  the  beginning  of  squeezing 
(Fig.  2.13),  the  radial  velocity  in  the  partially-filled 
case  is  somewhat  larger  than  in  the  fully-filled  case, 
especially  near  the  centerplane  at  the  edge  of  the 
disks(z*0,  r=R) .  This  can  be  understood,  considering  that 
the  fully-filled  case  experiences  some  resistance  from  the 
fluid  outside  the  edge  of  the  disks.  This  diffi^rence 
becomes  smaller  and  smaller  as  squeezing  goes  on  (see 
Fig.  2.14  2.16),  since  the  partially-filled  case  begins  to 
build  up  a  “bulge"  of  fluid  near  the  front,  which  acts  like 
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Read  data  (  X  ( I )  ,  FORCr.0 ,  ^  ,  etc  .  ) 

f  no 


Time  zero  ? 


Finite  element  solution  (V(I),  p(I)) 
Calculation  of  force  (FORCE) 


Linear  adjustment  of  solution 
V ( I ) =V ( I ) * ( FORCEO/FORCE ) 
p (I ) =p ( I ) * ( FORCED / FORCE ) 


Print  solution 


Next  time  step  ? 


'STOP 


^old 

J^ld(l)=V(l) 


Read  previous 
solution  (V(l),p(I)) 


Enter  At 


Move  boundary  nodes  (predictor  part) 
2new=5cad+Xold* 


Determine  new  contact  point 


Finite  element  solution  Pnew^^^^ 

Calculation  of  force  (FORCE) 

i 

Linear  adjustment  of  solution 

i 

Move  boundary  nodes  (corrector  part) 
Sna«=2old-"i(&ld  +X,ew  > 


Determine  new  contact  point 


Fig.  2.11  Numerical  algorithm  to  solve  the 
partially-filled  comprressive  flow  of  Newto 
nian  fluids. 
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R  =  5.042 
H  =  0.901 


R  =  5.517 
H  =  0.701 

Fig.  2.12  Movement  of  the  boundary  in  the  partially- 
filled  compressive  flow,  Ro/Ho=5* 


Fig.  2.13 
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1000. 
5.  042 
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vr/i 

Fig.  2.14  The  radial  velocity  profiles  of  the  partially- 
filled  compressive  flow  compaired  to  the  fully-filled  case 


Fig.  2.16  The  radial  velocity  profiles  of  the  partially- 
filled  compressive  flow  compared  to  the  fully-filled  case 
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the  fluid  outside  the  disks  in  the  full-filled  case. 
Pig.  2.17  shows  the  pathlines  of  some  material  points  as  the 
squeezing  goes  on.  The  material  points  near  the  front  are 
moving  upwards,  which  represents  the  same  phenomenon  as 
observed  in  the  injection  molding  process,  known  as  the 
"fountain  effect".  The  term  "fountain  effect"  was  coined 
and  discussed  by  Rose(1961)  and  is  important  in  determining 
the  quality  and  morphology  of  the  surface  of  the  molded 
article . 


The  movement  of  the  boundary  and  the  comparisons  of 
the  radial  velocity  profiles  are  shown  in  figures  2.13 
through  2.21  for  a  large  value  of  Ro/Ho(=15).  We  see 
similar  phenomena  as  for  Ro/Ho=5,  but  the  effect  of  the 
fluid  outside  the  disks  is  rather  small  in  this  case  and  it 
becomes  much  smaller  as  the  squeezing  goes  on  ( see 
Fig.  2.21).  We  can  therefore  conclude  that  as  long  as  Ro/H© 
is  large  enough,  the  flow  patterns  in  the  partially-filled 
compressive  flow  are  essentially  the  same  as  those  one 
observes  in  the  fully— filled  case  except  near  the  front: 
here  we  observe  fountain  flow  phenomenon  in  the 
partially-filled  case. 
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1^ 


Fig.  2.17  The  pathlines  of  material  points  in 
the  partially-filled  compressive  flow  of 
Newtonian  fluid  ;  the  fountain  effect  is  seen 
-  at  the  moving  front. 


Fig.  2.20  The  radial  velocity  profiles  of  the  partially- 
filled  compressive  flow  compared  to  the  fully-filled  case 


CHAPTER  3 


LUBRICATED  COMPRESSIVE  FLOW  OF 
VISCOELASTIC  FLUIDS 

3 . 1  Problem  formulation 

Let  us  consider  that  the  material  is  compressed 
between  two  circular  disks  under  a  constant  load  and  assume 
that  there  exist  thin  lubricant  or  low  viscosity  fluid 
layers  near  the  wall.  The  schematic  diagram  is  shown  in 
Fig.  3.1,  in  which  the  radius  of  the  disk  may  or  may  not 
vary  in  time.  If  the  viscosity  of  the  lubricant  is  given  in 
the  proper  range,  which  will  be  discussed  in  Section  3.5, 
the  flow  in  the  viscous  material  may  be  assumed  as  biaxial 
extensional  flow  with  negligible  shear.  The  velocity 
profile  in  the  cylindrical  coordinate  is  then  written  as 

Vy  =  I"  (3-1) 

V2  =  -26i,Cfe)Z  (3.2) 

\4  =  0  (3.3) 

in  which  is  the  biaxial  extensional  rate  and  varies  in 

time.  The  deformation  rate  tensor  d  is  given  by 
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z 


0  RltT 

(a) 


z 


0  R=Ro 


(b) 

Fig.  3.1  Schematic  diagram  of  the 
lubricated  compressive  flow. 

(a)  R  varies  in  time 

(b)  R=Ro  (constant) 
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where 


(s:-p  +  Tff)  is  the  radial  component  of  the  total 


47 


stress  tensor  ff;.  Evaluating  the  integration  in  (3.9), 

<iz 

=  -P.  +  7rf]  Jz 

=  +6b^)  “  --ZpoH  f  2  7rrH  =•  0 

Thus,  we  obtain  the  equation  for  po(t). 


p.(t)  =  +  V 


Substituting  into  (3.8)  then  gives 


p(r.z,'t)  =  ip(R‘-i-*)(#  +  4")-^P(«-3z“)(^-i4“)  +  7;r<3-ll 


The  total  force  exerted  by  the  fluid  on  the  upper 
dish  is  calculated  through  the  following  integration; 

F  =  f  f-Oi  ZA,rdr 

”  /  ^  f  “  "<zz  ^2  =  iH  (3.12) 

•'o 

Substituting  equation  (3.11)  into  (3.12)  and  performing  the 
integration,  one  obtains 

F=  +  +  |■7£/’«V(^-2e?j  +  /^R^(7^^-7zz)  (3.13) 


One  of  the  easiest  ways  to  apply  the  force  on  the 
upper  dish  is  to  place  a  mass  on  it.  In  that  case,  the 
force  is  given  by 


F  =  mi  3-a) 


(3.14) 


where  m  =  the  total  mass  on  and  in  the  disk 
g  =  the  gravitational  acceleration 
a  =  the  acceleration  of  the  disk  i  -  -2. 


Equivalently, 


m(  3  -4H-^  +  8h4'^) 


(3.15) 


(Vz)z  =  H  = 


(3.16) 


a  =  -2 


±/dH 


m 


=  -  8H€j,^ 


(3.17) 


Now,  equating  equations  (3.13)  and  (3.15),  we  have 

(3-+h#*8H67)  =  + 


+  ~  '^zz.  ) 


(3.18) 


mid  +  8 Hii!")-  x7l/>r^4^  +  -  7at( ) 


K*[/  +  />x{xy  ■  rcR 
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Frequently,  equations  (3.18)  or  (3.19)  can  be  simplified  as 
follows : 


(a)  When  R>>H,  the  term  including  H  on  the 
righthand  side  of  (3.13)  can  be  neglected  compared  to  the 
first  term  to  give 


c/4 

d-t 


m(3+8Hik^)-  y».- ) 


(3.2o; 


(b)  When  the  compression  is  not  very  fast,  the 
inertia  of  the  fluid  and  the  load  can  be  neglected,  in  which 
case  equation  (3.18)  is  simplified  to 

IflJ  ~  TiR  (  Tf-r  ~  )  (3.21) 


When  the  radius  of  disk  varies  in  time  (see 
Fig.  3.1(a)),  the  expression  for  the  radius  change  is 
obtained  from  mass  conservation;  that  is. 


or 


or 


=  constant 


=  0 


(3.22) 


Here  aquation  (3.16)  has  been  used. 


The  next  step  is  to  solve  equation  (3.19)  or  the 
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simplified  form  ( ( 3 . 20 ) , ( 3 . 21 ) )  together  with  (3.16), 
(3.22),  and  appropriate  constitutive  equations  which  relate 
Trr  Tzz  to  the  deformation.  In  the  next  section,  we 

will  consider  the  linear  viscoelastic  case,  since  in  the 
limit  of  small  deformations  the  response  of  polymeric 
materials  may  be  considered  linear.  The  non-linear  case 
will  be  considered  in  Section  3.3. 

3 . 2  A  linearized  viscoelastic  case 

For  sufficiently  small  values  of  the  deformation, 
the  mechanical  behavior  of  polymeric  materials  is  entirely 
described  by  the  constitutive  equation  of  linear 
viscoelasticity . 

3.2.1  Constitutive  equation 

A  particular  case  of  the  constitutive  equation  of 
linearized  viscoelasticity  is  given  by 

2  +  (3.23) 

where  Ai  stands  for  the  relaxation  time  and  stands  for 

the  retardation  time.  Often  Aa  is  set  to  be  zero.  Note 
that  as  Aa  approaches  A| ,  equation  (3.23)  approaches  the 
Newtonian  equation.* 


*  When  Ai=Ai  ,  aquation  (3.23)  is  simply  the  sum  of  the 
Newtonian  constitutive  equation  and  the  first  time 
derivative  multiplied  by  Aj. 
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In  the  present  problem,  the  radial  and 
components  of  the  constitutive  equation  are 

7,, +  A,^=a^(4  +  A,#) 

T;»  +  =  -4)  (  4  + 

from  which  it  follows  that 

"^zz.  ~  ~  ^  'Ti-ir 

and 

=  exp(- ^j'2- jj- 1  ( 4  +  Az  exp  ( -^ ]  dt 

O 

3.2.2  Analytical  solution 

I 

Let  us  consider  the  case  ( b)  in  Fig.  3.1,  in 
the  radius  is  constant,  and  that  R=Ro>>H.  Further 
that  the  compression  is  not  very  fast,  so  that 
that  the  inertia  of  the  load  is  small  compared 
inertia  of  fluid  (  3.20).  Then  equation 

or  (3.20)  becomes 

or 

dei,  _  A- m3  _  12. 

dt  “  7Cf  izj-  PR^- 

Here  equation  (3.26)  has  been  used.  Substituting 
into  (3.28)  then  gives 


axial 

(3.24) 

(3.25) 

(3.26) 

(3.27) 

which 
assume 
-  ,  and 
to  the 
(3.13) 

(3.28) 

(3.27) 


dt  “  /’Ro'^A,  ^  t) 


Differentiating  this  with  respect  to  time  gives 


■^p(i)#+3r-''Kxj^  - 


Tp^x: 


■2.4  ^Ai  \  c/^fc 


which  is  a  second  order  dynamical  equation 


The  solution  of  equation  (3.29)  is  given  by 


J^i  g.  c,exp(D,tJ  +  C:,exf>(i>^-t) 

^  “  ivWn  D,%D^ 

+«p(D,t).  (Cjt  +C,) 

'  ivhen 


(3.30) 


where  D;  and  are  the  roots  of  the  characteristic  equation 


of  (3.29),  that  is, 


/>A,P^  +  />(  I  +  )  0  +  =  0  (3.31) 


Cf  t  Cz,  C3 ,  C4.  are  constants  to  be  determined  by  the 

following  initial  conditions; 
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The  characteristic  feature  of  the  solution  (3.30) 
may  be  divided  into  the  following  three  cases,  depending 
upon  the  values  of  D,  and  • 

Case  I.  (I  -  >  0  (3.33) 

The  Dj  are  real,  negative,  unequal  and  given  by 

i  [-  0+ 

The  solution  decays  exponentially  to  a  constant  rate;  that 
is, 

•  (Y)  Cf 

^  C,^p(P,t)  +  C;,€xp(P,tJ 

where 

C,  = 

c,  =  — ! _ 1 0 ..  4-  -AulL.] 

Note  that  the  solution  approaches  the  Newtonian  solution  as 
Ai  approaches  to  Xi  >  long  times  for  any  A|  and  A^. 
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Newtonian  solution  is  given  by 

m  3 


Case  II.  ( 


_  _  9gn/  _ 


Pi^i- 


-  0 


(3.35) 


P I  and  are  equal  and  negative. 


D,  =  Di  =  - 


-^Aj 


and 


where 


■  (4*+  C,) 


(3.36) 


Ca  — 


AmS 


npizj-  5^7 /2o-^ 

m3 


md 


^4  = 

That  is,  the  fluid  is  gradually  accelerated  to  a  steady 
state  deformation  rate. 

Case  III.  (  I  +  - ^^IAl 


pRP-  •  PRP- 

Di  and  Dz  complex  conjugates,  given  by 


<  0 


(3.37) 


and  the  solution  is  an  oscillatory  one  in  which  we  are 
particularly  interested.  The  following  case  is  of 
particular  interest: 


(3.38) 
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and,  from  (3.16) 


Each  term  in  equation  (3.40)  represents  different  features 
of  the  lubricated  compressive  flow  of  linear  viscoelastic 
materials  under  the  condition  (3.38)  .  The  first  term  (A.) 


represents  Newtonian  viscous  damping,  the  second  (B)  is 
related  to  the  initial  amplitude,  and  the  third  (C) 
represents  the  damping  of  the  oscillatory  term  (D),  whose 

ITT' 

oscillation  period  is  7lRo\j  • 


This  oscillatory  motion  originates  from  the 

9  V 

coexistence  of  the  unsteady  inertia  term  and  the 

ah’ 

unsteady  elastic  relaxation  term  (A  •  These  two  time 
derivatives  interact  to  generate  the  second  order  time 
derivative  in  equation  (3.29). 


The  condition  (3.37)  is  rather  conservative  for  the 
oscillatory  motion,  since  the  solution  (3.40)  will  behave  as 
the  non-oscillatory  one  if  the  damping  term  (C)  is  fast 
enough  to  make  the  product  (C)*(D)  very  small  in  a  time 
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which  is  smaller  than  the  oscillation  period.  From  this 
fact,  we  can  develop  a  more  useful  necessary  condition  for 
the  oscillatory  motion:  that  is, 

time  scale  of  damping  >  oscillation  period 


It  is  useful  to  rewrite  the  condition  (3.42)  and 
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(3.43)  in  terms  of  the  dimensionless  groups  defined  below: 


= 

oscillation 

(3.44) 

oT  = 

no  oscillation 

(3.45) 

Hyp 

(Reynolds 

No. ) 

Re  = 

7 

^  V  Aj 

=  H 

(Deborah 

No. ) 

V  Xz 

Rd=  ^ 

(Retardation  No.) 

Note  that  or  depends  on  the  material  properties 
and  the  geometry  only,  but  not  on  the  squeezing  speed,  even 
though  each  group  does  depend  on  the  squeezing  speed.  Also 
note  that  the  dimensionless  group  ^®/Re  similar  to  the 
first  elasticity  number  given  by  Astarita  and 
Marrucci(1974) ,  which  has  been  defined  by  the  ratio  of  the 
Weissenberg  and  Reynolds  numbers  in  quasi-viscometric  flows. 
This  elasticity  number  was  employed  by  Denn  and 
Porteous( 1971 )  to  identify  conditions  under  which  elasticity 
can  be  expected  to  be  important  in  viscoelastic  fluid  flow. 
Also  it  is  seen  in  the  paper  by  Tordella (1953 )  who  used  this 
number  to  predict  the  onset  of  melt  fracture,  and  in  the 
paper  by  Soger (1977)  who  tried  to  correlate  pressure  losses 
due  to  the  elasticity  in  the  capillary  rheometer  to  the 


58 


elasticity  number. 

The  equation  (3.29)  and  the  conditions  (3.42-3.45) 
can  be  laodified  to  include  the  load  inertia  terra,  which  has 
been  neglected  in  the  above  by  assuming  First, 
we  define  several  different  densities. 


P  ^  Pm 


(3.46) 


m  H 
71 


(3.47) 


yo  is  the  fluid  density,  and  is  a  pseudo-fluid  density 
corresponding  to  the  mass  of  the  load  (m)  in  the  lubricated 
compressive  flow;  is  the  apparent  density  which 
represents  both  effects  of  the  fluid  inertia  and  the  load 
inertia.  If  we  use  instead  of  using  p  in  the  equation 
(3.29)  to  include  both  inertia  effects,  we  have 


(3.48) 


Here  the  constant  radius  condition  also  has  been  relaxed  to 
include  the  situation  of  Fig.  3.1(a).  Then  the  conditions 
(3.42-45)  become 

>  ( I  +  '  oscillation  (3.49) 

<  (i  P  :  no  oscillation  (3.50) 

t  K  r  fS 


or 
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;  oscillation  (3.51) 


:  no  oscillation  (3.52) 


where 


Ra  = 

1 


(modified  Reynolds  No.) 

I 

The  oscillation  period  is  also  modified  to  JcRy  . 


At  this  point  we  may  note  several  interesting 
aspects  of  the  lubricated  compressive  flow  of  linear 
viscoelastic  materials  ; 


(1)  The  compressive  motion  under  a  constant  amount 
of  load  may  or  may  not  be  oscillatory,  depending  upon  the 
conditions  (the  material  properties,  the  amount  of  the  load, 
and  the  geometry)  as  predicted  by  equations  ( 3 . 49 )- ( 3 . 52 ) . 


(2)  When  oscillation  occurs,  it  is  due  to  the 
combined  affects  of  the  inertia,  the  elasticity,  and  the 
viscosity.  The  oscillation  period  depends  only  on  geometry, 
inertia,  and  the  elastic  modulus. 


(3)  The  retardation  time  may  play  an  important  role 
in  the  damping  of  the  oscillatory  motion  (see  term  (C)  in 
equation  (3.40)). 

(4)  As  Ai  approaches  A|  ,  the  response  of  linear 
viscoelastic  materials  approaches  the  Newtonian  one. 
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In  the  next  section,  numerical  solutions  will  be 
considered  with  no  assumptions  made,  and  this  will  further 
support  and  refine  the  conclusions  noted. 


3-2.3  Numerical  solution 

Equation  (3.19),  together  with  (3.16),  (3.22),  and 
constitutive  equations  (3.24)  and  (3.25),  can  be  solved 
simultaneously  by  the  Runge-Kutta  method  or  the  method  of 
Gear(1970) .  Many  different  combinations  of  material 
properties  and  geometrical  factors  have  been  considered  to 
investigate  the  characteristic  features  of  this  flow  in  the 
case  of  Pig.  3.1(a),  in  which  the  radius  varies  to  keep  the 
total  mass  of  the  fluid  constant. 


In  Fig.  3.2,  curves  of  H/Hq  vs.  time  are  shown  at 
various  values  of  relaxation  times  (  A|  =  0.003,  0.03,  0.1, 
0.3)  and  zero  retardation  time,  and  compared  to  the 

corresponding  Newtonian  curve  (  A,  =  0.).  When  to . is 

less  than  the  unity,  the  curve  does  not  oscillate  (see  curve 


4),  as  expected  from  the  condition  (3.50).  As  /O 


lA, 


increases,  the  oscillatory  motion  begins  and  becomes  severe 
as  the  relaxation  time,  and  hence  the  value  of  the  group 
to  ,  increases  (see  curves  1  to  3  ).  The  oscillation 
period  is  about  the  same  as  the  calculated  value  from 
aquation  (3.40),  if  we  consider  the  change  of  radius  in 
these  calculations  [equation  (3.40)  was  derived  on  the  basis 


00  * 


flow  of  linear  viscoelastic  materials  ;  the  effect  of 
the  relaxation  time. 


the  retardation 
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of  a  constant  radius]  .  Note  that  itiost  of  the  viscoelastic 
curves  remain  below  the  corresponding  Newtonian  curve, 
though  the  difference  is  minor  at  large  times. 

In  Pig.  3.3,  the  effect  of  the  retardation  time  has 
been  illustrated.  All  the  conditions  are  the  same  as 
before,  except  that  X,  is  kept  constant  (at  O.l)  in  all 
these  calculations,  and  Xa  varies  (  Xi  =  0.,  0.01,  0.05, 
0.1).  We  can  recognize  that  a  small  increase  of  the 
retardation  time  reduces  the  oscillation  amplitude  greatly, 
and  the  response  approaches  the  Newtonian  curve  as  Aa 
approaches  Xj.  All  curves  remain  below  the  corresponding 
Newtonian  curve  (X|=Xa=0.). 

3 . 3  Non-linear  viscoelastic  cases 

3.3.1  Constitutive  equations 

Since  the  linear  model  does  not  predict  some  of  the 
important  non-Newtonian  behavior  of  polymeric  materials, 
such  as  shear  thinning,  normal  stresses  in  simple  shear 
flows,  and  high  extensional  viscosity  in  extensional  flows, 
various  kinds  of  non-linear  models  have  been  proposed.  We 
snail  now  consider  some  of  the  currently  popular  and 
promising  models. 

One  of  the  most  successful  models  is  the  one 
proposed  by  Phan-Thien  and  Tanner (1977)  and 
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Phan-Thian( 1978 ) ,  which  has  the  following  form  with  a  single 
relaxation  time; 

T  1]  +  A[(l-i;I)T  +  7  ]  =  i  (3.53) 

Hera,  7  and  7  are  the  contravariant  and  •  covariant 
convected  derivatives  given  by  Oldroyd ( 1950 ) ,  respectively, 
and  denoted  by 

r  =  -7V'r-T 

A  n7  T 

X=  -fe+vvT-hivy 

Dt  ~  ~  ~  ~ 

This  equation  has  four  material  parameters,  X  »  5  »  and 

£  .  %  is  the  zero  shear  viscosity  and  A  is  the  relaxation 

time.  The  parameter  3  t)e  determined  from  normal  stress 

data  in  simple  shear  flow,  since  is  simply  the  ratio  of 

the  second  normal  stress  difference  to  the  first  normal 
stress  difference,  or  from  the  deviation  of  the  viscosity 
from  8  is  related  to  the  limiting  value  of  extensional 

viscosity  by  ^  =  const  -  (^/s) 

7g-*90 

When  8-  vanishes,  (i.e.,  when  the  extensional 
viscosity  has  no  limit) ,  the  Phan-Thien-Tanner  model  reduces 
to  the  one  by  Johnson  and  Segalman( 1977 ) ,  represented  by 


y  +  Afd-is)!  *  ssil  = 


(3.54) 
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By  adding  a  retardation  term  to  the  righthand  side  of 
(3.54),  a  modification  of  the  Johnson-Segalman  model  is 
obtained  as 

* 

where  A^,  is  the  retardation  time  and  M  and  a  are  defined 
in  the  same  way  as  “f  and  7  -  The  existence  of  the 
retardation  terra  can  be  explained  if  one  considers  a  polymer 
solution  made  up  of  a  Newtonian  solvent  and  a  polymeric 
solute  which  is  described  by  (3.54).  The  total  stress  in 
the  solution  would  be  the  sum  of  the  polymeric  contribution 
and  Newtonian  solvent  contribution,  that  is. 


(3.57) 


(3.53) 

Here  the  subscript  p  stands  for  the  polymer  and  sv  for  the 
solvent . 

Differentiating  equation  (3.53) 

Xs*  -  ^  40(1  Tsv  * 

or 

(|'£5?Ai2sv  = 


(3.59) 
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I  r  ^  ^  1  o 

*  it  X  2J  is  often  denoted  by  X  is  known  as 

corotational  Jaumann  derivative  (see  Zaremba( 1903 )  and 
Fromm ( 1947 ) ) . 
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where  is  the  second  invariant  of  the  deformation  rate 
tensor  and  G  is  the  shear  modulus.  This  model  is  useful  in 
situations  in  which  the  non— Newtonian  viscosity  plays  an 
important  role . 

The  model  proposed  by  Acierno  et.  al. {1976a),  often 
called  the  "structural  model",  is  also  very  promising.  The 
modal,  with  a  single  relaxation  time,  is  represented  by 


'1  +  2)  = 


(3.67) 


^  X  ,  A  =  Ao 

^  (a=  0.25 -0.4} 

E  = 

where  is  the  contravariant  convected  time  derivative, 

Oc 

and  the  scalar  dimensionless  quantity  x(  £1.)  is  regarded  as 
a  structural  variable  which  describes  how  far  the  existing 
structure  is  from  equilibrium.  G©  and  A©  the 

equilibrium  values  of  G  and  A  (when  x=l.). 


For  comparison,  we  summarize  the  responses  of  each 
model  fluid  in  simple  shear  flow  in  Table  3.1. 


The  covariant  convected  Maxwell  model  predicts  too 
large  a  value  of  Vii  -  laa  in  comparison  to  what  is  found  in 
most  polymeric  materials.  The  corotational  model  exhibits  a 
shear  viscosity  which  depends  upon  the  shear  rate  too 
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^iX 

“  ^22. 

Stress 

overshoot 

Phan-Thien- 

Tanner 

No  exploit  expression,  but 
very  close  to  Johnson-Segalman 

yes 

Johnson- 

Segalman 

Xf 

zXXlr^ 

-lU?" 

yes 

1  +A^fif2-5) 

Johnson-Segalman 
with  retard,  time 

2)o(K-X2)i^ 

2r 

yes 

I  +x’li^s(2-i) 

Contravariant 
convected  Maxwell 

0 

no 

Covariant 
convected  Maxwell 

-z^Xr"" 

no 

Corotational 

Maxwell 

')o  y 

-\xf^ 

yes 

/+A^P 

b&h 

1  +  x^f^ 

Vfhite- 

Metzner 

5 

0 

no 

Marrucci  ^ 

structural 

0 

yes 

I  X  A  * 

*  Structural  variable  x  satisfies  - Jr?~  “  AAo^' 

where  ^1=  0. 25 -^0 . 4.  ^ 


Table  3.1  Comparison  of  non-linear  models  in 
simple  shear  flow. 
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strongly;  it  also  pradicts  spurious  oscillations  after 
stress  overshoot  and  much  too  large  a  value  of  Tli  ”  • 

Therefore,  we  shall  disregard  these  two  models  in  our 
numerical  calculations. 

The  contravar iant  convected  Maxwell  model  and  its 
modification,  the  Whits— Metzner  model,  predict  zero 
and  no  stress  overshoot  during  startup.  The 

Phan-Thien-Tanner  modal  and  its  simplified  form,  the 

Johnson-Segalman  model,  predict  non-zero  values  of 
and  also  a  stress  overshoot.  The  structural  model  predicts 
zero  but  it  does  predict  a  stress  overshoot. 

3.3.2  Numerical  calculations 

Lubricated  compressive  flows  of  the  viscoelastic 
model  fluids  mentioned  above,  including  the  contravariant 
convected  Maxwell  model,  the  White— Metzner  model,  the 
Johnson-Segalman  model  {with  and  without  retardation  term) , 
and  the  structural  model,  are  solved  numerically  in  this 
section.  Since  the  aquations  to  be  solved  in  each  case, 
which  are  given  in  appendix  B,  consist  of  first  order 
ordinary  differential  equations,  the  Runge-Kutta  or  Gear 
method  can  be  used.  Here,  we  consider  the  situation  given 
in  Fig.  3.1(a),  and  the  initial  conditions  are  given  by 


The  results  are  shown  in  Figs.  3.4  through  3.10,  in 
which  H/Hq  is  plotted  against  time  in  each  case.  The  curves 
for  the  contravariant  convected  Maxwell  fluid  are  given  in 
Fig.  3.4  and  3.5,  those  for  the  White-Met zner  model  in 
Fig.  3.6,  the  Johnson-Segalman  model  in  Figs.  3. 7-3. 9,  and 
the  structural  model  in  Fig.  3.10. 


When  _1A_  is  small  enough,  the  oscillation  does 

^  ^  7  A 

not  occur,  as  shown  in  Fig.  3.4.  But  if  is 

small,  oscillation  does  occur.  When  the  squeezing  rate  is 


ilow  or  moderate,  the  oscillatory  curves  stay  below  the 


corresponding  inelastic  ones  in  all  cases  (see  Fig.  3.5, 
3.6.,  3.7,  3.10).  In  fast  squeezing,  each  model  behavior  is 


distinct . 


The  Maxwell  model  results  in  oscillation  and  stays 
below  the  Newtonian  curve  (Fig.  3.5).  In  contrast,  the 
White-Metzner  model  does  not  give  oscillation  at  all  for  the 
power— law  index  of  0.5  and  stays  below  the  corresponding 
power-law  curve  (Fig.  3.6).  The  absence  of  oscillation  in 
the  White-Metzner  model  in  fast  squeezing  seems  to  be  due  to 


the  large  decrease  of  the  viscosity  and  the  relaxation  time 


at  high  deformation  rates.  These  effects,  plus  the 
increases  in  R  as  time  progresses,  combine  to  make 

r  ^ 

very  small,  under  which  condition  no  oscillation  would  be 


expected . 


71 


The  faster  squeezing  of  the  Maxwell  model  may  look 
unusual  to  those  who  consider  that  the  Maxwell  model 
predicts  very  large  elongational  viscosities.  But,  if  one 
remembers  correctly,  the  large  elongational  viscosities  of 
Maxwell  fluids  are  obtained  at  large  values  of  after  a 
long  time.  In  other  words,  in  a  short  time  or  at  small 

values  of  4A  the  elongational  viscosity  is  even  smaller 

• 

than  in  the  Newtonian  case  (5|).  In  the  squeezing  flow,  €1, 
is  large  initially,  but  this  large  4  maintained 
for  a  long  time  since  the  film  thickness  is  getting  thinner 
very  quickly.  Thus,  it  is  very  difficult  to  build  up  a 
large  elongational  viscosity  in  this  geometry.  This  is  why 
the  Maxwell  fluid  is  squeezed  faster  than  the  Newtonian 
fluid  even  with  the  high  initial  value  of  4A.  But,  under 
the  extraordinarily  high  loading  conditions,  the  Maxwell 
curves  are  crossing  over  the  corresponding  Newtonian  ones  at 
very  small  values  of  H/Hq  (see  \ppendix  G  for  the  detail). 

Even  though  the  absolute  differences  between  (H/Hq  )  Maxwell 

-3 

and  (H/Ho)M«wton;an  small  (in  the  order  of  10  or 
less) ,  the  relative  differences  between  them  are  around  2  in 
the  cases  of  Appendix  G.  This  fact  may  be  important  in  the 
field  of  lubrication  technology. 

The  rapid  squeezing  of  the  Johnson-Segalman  model 
(Fig.  3. 7-3.9)  and  the  structural  model  (Fig.  3.10)  are 
particularly  interesting,  because  those  curves  move  above 


fluid  :  no  oscillation  at  small 


00* 


flow  of  contravariant  convected  Maxwell  model,  compared 
to  Newtonian  fluid. 


flow  of  Johnson-Segalman  model 
responding  shear  thinning  case 


00* 


flow  of  Johnson-Segalman  model  with  non-zero  retarda¬ 
tion  time,  compared  to  the  corresponding  inelastic  case 


INELASTIC 
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the  corresponding  inelastic  ones  after  some  amount  of 
initial  faster  squeezing.  The  effect  is  more  prominent  in 
the  Johnson-Segalman  model.  The  slower  squeezing  of  these 
models  must  be  due  to  the  specific  features  which  are  not 
found  in  other  models:  that  is,  the  predictions  of  the 
stress  overshoot  phenomenon  and  non-zero  -  Tag  in  simple 
shear  flow  of  the  Johnson-Segalman  model,  and  the  stress 
overshoot  in  the  structural  model.  Figs.  3.3  and  3.9  show 
the  effect  of  non-zero  retardation  time  in  the 
Johnson-Segalman  model,  which  damps  the  oscillation. 

3.4  Finite  element  simulation  (Maxwell) 

The  main  purpose  of  this  section  is  to  develop  and 
test  a  finite  element  numerical  scheme,  to  be  used  in 
solving  the  unlubricated  problem  in  chapter  4.  In 
particular,  we  will  consider  the  contravariant  convected 
Maxwell  model.  The  problem  has  been  already  solved  in  the 
lubricated  case  in  Section  3.3.2.  Thus,  the  numerical 
scheme  is  tested  by  solving  that  same  problem  and  comparing 
the  result  to  the  one  obtained  in  Section  3.3.2. 

The  problem  becomes  complicated  since  we  wish  to 
include  the  unsteady  inertia  and  the  unsteady  elasticity 
terms  in  the  equations  of  motion  and  the  “Constitutive 
equations,  respectively.  Numerically,  it  becomes  more 
complicated  because  the  domain  changes  as  time  goes  on  and 
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unsteady  transient  equations  have  to  be  solved  between  time 
steps.  Thus,  we  will  pose  the  problem  in  a  more  natural 
way,  following  the  pathlines  of  each  material  point  in 
Lagrangian  coordinates.  In  other  words,  each  nodal  point  in 
the  initial  grid  moves  along  its  own  pathline  as  time  goes 
on,  and  we  use  material  time  derivatives  instead  of  using 
the  partial  time  derivatives  in  the  equations. 

The  equations  to  be  solved  are 

7‘V  =  0  (3.68) 

dV 

=  -VP  +  7-2  (3.69) 

2  +  -vyT-2vy’^]  = +  (3.7o) 

Initial  conditions  have  to  be  specified  (since  we  have  time 
derivatives) .  These  initial  conditions  are  given  by 


r  = 

y:  =  0  at 

t  =  a 

(3.71) 

The 

time 

derivatives  are 

treated  by  a 

finite 

difference 

scheme 

in  the  time 

coordinate.  We 

use  the 

implicit  three  point  recurrence  scheme  with  variable  time 
steps,  which  requires  at  least  two  previous  solutions.  To 
solve  for  the  solution  at  t=tn^.,  based  on  two  previous 
solutions  (at  tn  and  tn-i  )  ,  the  material  time  derivative  is 
discretized  (on  the  time  axis)  by 
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- L - [  4t„-,  (24tn  +  ^'tn-i)A  nfl 

AZn  (  ^tn +^^*1-.!  )  u 

-  (Atfi  +  An  +(^tn)  >An-il^ 2  72) 

where  At„.,  =  tn  -  tn-i 

jfl'fcn  =  "fen+i  “  "fcn 

A[  =  A  i  —  ii 


All  other  terms  in  the  equations  involve  the  unknown 
variables  at  t=tn+i  •  The  discretized  (in  time)  forms  of 
equations  ( 3 . 68 ) - ( 3 . 70 )  are  given  by 


^  2-. 

I.„  +  =  1  [«."  +  ”""]  (3.75) 

The  whole  system  has  to  be  solved  simultaneously  on  the 
given  domain.  The  Newton-Raphson  method  is  used  to  treat 
the  nonlinear  terms,  and  the  predictor-corrector  method  is 
used  to  improve  the  shape  of  the  grid  in  each  time  step. 


The  mixed  finite  element  method,  which  was  first 
introduced  by  Kawahara  and  Takeuchi( 1977 )  and  further 
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studied  by  Crochet  and  Bezy(l979),  and  Crochet  and 
Keunings ( 1980) ,  has  been  adopted  to  solve  the  system 
(3 . 73 )- ( 3 . 75 ) .  This  method  differs  from  the  displacement 
method  used  by  Chang  et.  al.( 1979a, b)  in  terms  of  the 
unknown  fields.  In  the  mixed  method,  the  pressure,  the 
velocity  components,  and  the  stress  components  are  the 
unknown  fields,  while  in  the  displacement  method,  the 
pressure  and  the  velocity  components  are  the  unknown  fields 
and  the  stress  components  are  calculated  by  means  of  an 
iterative  technique.  In  particular,  we  have  chosen  the 
triangular  elements  and  the  same  shape  functions  as  those  of 
Crochet  and  Keunings;  that  is,  the  linear  shape  function  in 
the  pressure  and  the  quadratic  shape  function  in  the 
velocity  and  the  stress  components  (see  Fig.  3.11).  The 
Galerkin  finite  element  formulation  of  system  ( 3 . 59 ) - ( 3 . 61 ) 
is  straightforward  and  is  given  in  Appendix  D.  The 
resulting  simultaneous  linear  algebraic  system  is  solved  by 
the  frontal  elimination  technique  proposed  by  Irons(1970)  to 
reduce  the  use  of  the  central  core  memory.  The  numerical 
algorithm  to  solve  the  lubricated  compressive  flow  of  the 
contravariant  convected  Maxwell  fluid  under  a  constant  load 
is  given  in  Fig.  3.12. 


The  initial  grid  used  in 

the 

calculation 

and 

the 

deformed  grids  at  later  times 

are 

shown  in  Fig . 

3.13, 

in 

which  one  can  easily  see  that 

the 

deformation 

of 

the 
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Read  previous  solutions  (  Yn*2n>  Pn  ) 

I  (  Yn-l,rn-l.  Pn-1  ) 


- >  Enter  At 

i 

Move  nodal  points  (  X^+l  ) 


Finite  element  solution  (  Vn+l»2n+l»  Pn+1  )< — i 


Calculation  of  force 

I  Adjust  top  disk 

Convergence  test  ? - ^velocity 

I  yes 

Print  solution 

^  no 

Next  time  step  ? - ^STOP 

I  yes 

_ Yn-l~Xn»  2’n-l~2n»  Pn-l“Pn 

XntYn+1 «  Jn  =  Zn+l »  Pn“Pn+l 


Fig.  3.12  Numerical  algorithm  to  solve 
the  lubricated  compressive  flow  of  Maxwell 
fluid. 


H/Ho  =1.0 
Time  =  0.0 


H/Ho  =  0.7641 
Time  =  0.075 


H/Ho  =  0.8198 
Time  =  0.120 


H/Ho  =  0.6938 
Time  =  0 . 230 


Fig.  3.13  The  initial  grid  and  the  deformed 
grids  (at  later  times)  in  the  lubricated  com¬ 
pressive  flow  of  contravariant  convected 
Maxwell  fluid. 
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material  is  indeed  biaxial  extension.  Small  numbers  of 
elements  have  been  used  in  this  calculation  because  the 
stress  field  is  uniform  throughout  the  domain,  and  there  is 
no  stress  singularity  in  this  problem.  Pig.  3.14  shows  the 
oscillatory  finite  element  result  (H/H^  vs.  time  curve) 
compared  to  the  semianalytical  solution  obtained  in  Section 
3.3.2.  The  agreement  between  them  is  very  good.  Now,  we 
can  proceed  to  the  next,  more  complicated  problem,  which  is 
the  unlubricated  compressive  flow  of  the  contravariant 
convected  Maxwell  fluid. 

3.5  Viscosity  range  of  the  lubricant  to  produce  the  biaxial 
extensional  flow 

If  the  viscosity  of  the  lubricant  is  too  low,  the 
lubricant  layer  will  be  expelled  quickly  during  the  early 
time  of  the  squeezing.  On  the  other  hand,  if  the  viscosity 
is  too  high,  then  the  flow  in  the  viscous  testing  fluid  will 
no  longer  be  extensional  flow.  Therefore,  to  generate  a 
lubricated  compressive  flow,  or  a  biaxial  extensional  flow, 
a  lubricant  with  the  proper  range  of  viscosity  has  to  be 
used . 


We  know  from  Section  2.2  that  the  less  viscous  fluid 
near  the  wall  is  expelled  preferentially  when  the 
extensional  stress  in  the  central  viscous  fluid  is  greater 
than  the  shear  stress  in  the  less  viscous  fluid  near  the 
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wall.  In  other  words,  in  order  to  ensure  that  the  lubricant 
layer  is  not  preferentially  expelled  during  the  early  time 
of  the  squeezing,  we  require  the  following  condition: 

extensional  stress  in  the  viscous  fluid 


shear  stress  in  the  lubricant 


<  1 


or,  for  Newtonian  fluids. 


4  ^  of  < 


(3.76) 


where  is  the  viscosity  of  the  lubricant,  ')x  is  the 
viscosity  of  the  central  viscous  material,  and  S  is  the 
thickness  of  the  lubricant  layer.  From  this  we  have 


(3.77) 


which  gives  us  the  lower  limit  of  the  lubricant  viscosity. 


If  the  condition  (3.77)  is  satisfied,  parallel 
squeezing  would  be  expected  and  hence  the  radial  velocity 
profile  is  given  by 


where 


(3.78) 


(3.79) 


See  Appendix  A  for  the  details  of  deriving  the  above 
velocity  profile. 
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In  order  to  have  nearly  extensional  flow  in  the 
central  viscous  fluid,  the  radial  velocity  at  the  interface 
between  two  fluids  (at  z=H-5)  should  be  very  close  to  the 
one  at  the  centerplane  (at  z=0).  Thus  we  can  set  a 
criterion  for  this  purpose  to  be 


=  -  >  0.95  (3.80) 

(Vr)z=o 


Therefore,  the  conditions  (3.77)  and  (3.81)  together 
determine  the  range  of  the  viscosity  of  the  lubricant  which 
should  be  used  to  produce  a  nearly  biaxial  extensional  flow 
in  the  central  viscous  fluid,  to  the  extent  that  the 
rheology  of  both  fluids  can  be  taken  to  be  Newtonian. 


CHAPTER  4 


ONLUBRICATED  COMPRESSIVE  FLOW 
OP  VISCOELASTIC  MATERIALS 

We  consider  the  unlubricated  compressive  flow,  in 
which  the  no-slip  boundary  condition  is  satisfied  by  the 
viscoelastic  materials.  First,  we  consider  a  linearized 
viscoelastic  case,  in  which  an  approximate  solution  can  be 
obtained,  which  gives  us  much  useful  information.  The 
contravariant  convected  Maxwell  fluid  will  then  be  examined, 
since  this  model  is  on  the  side  of  convenience  with  a  modest 
approach  toward  realism,  even  though  it  cannot  predict  the 
stress  overshoot  and  the  second  normal  stress  difference  in 
shear  flow.  In  this  case  a  numerical  solution  of  the 
partial  differential  equations  is  inevitable.  A  finite 
element  method,  which  has  been  tested  in  Section  3.3,  will 
be  used. 

4.1  A  linearized  viscoelastic  case 

4.1.1  Assumptions  and  governing  equations 

First  it  will  be  assumed  that  the  velocity  profiles 
are  the  same  as  those  of  a  Newtonian  fluid;  that  is. 
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Vz  = 


(4.1) 


Vj-  =  -^r  f'(z,i) 


(4.2) 


where 


(4.3) 


The  prime  in  equation  (4.2)  implies  partial  differentiation 
with  respect  to  z.  The  term  V  in  (4.3)  is  the  downward 
velocity  of  the  top  plate.  The  coordinate  system  used  here 
is  given  in  Fig.  4.1. 

We  will  also  assume  that  the  fluid  inertia  is  much 
less  important  than  the  load  inertia.  This  will  be  shown 
later  to  be  observed  in  most  of  the  experimental  conditions. 

From  equations  (4.1)  and  (4.2),  the  deformation  rate 
tensor  is  given  by 


j-ir  o 

^  =  -irf"  -if'  0 


(4.4) 


Using  the  constitutive  equation  (3.23),  the  non-zero 
components  of  the  stress  are: 


(4.5) 


:4.6) 


Fig.  4.1  The  domain  used  in  the  unlubri¬ 
cated  compressive  flow  of  viscoelastic 
materials . 
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■?iz  +A,^  =  i^(f'+A.|^) 
t 


(4.7) 


(4.3) 


from  which  each  component  of  the  extra  stress  tensor  7] 


obtained  as 


(4.9) 


7=  =  -2  7,. 


(4.10) 


^rz  =  -iT-'"exp(-^jJ  (4.11 


From  equation  (4.9)  ^  =  0  and  7^^-7e*=  0.  Therefore,  the 
equations  of  motion  become  (neglecting  fluid  inertia): 


(4.12) 


'0 


(4.13) 


(4.14) 


(4.15) 


Integrating  equations  (4.14)  and  (4.15),  the  pressure  is 
obtained  as  ^ 

p(r,2,t)  =-llr^exp(-±)  | 

f  +  poW(4.16) 

^  Q 
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To  determine  p^{t) ,  the  boundary  condition  at  the  free 
surface  (at  r=R)  is  applied  as 


(4.17) 


We  thus  have 


+2iexp(-j7)|*(V  +  AaiT)exp(^)<ft  (4.18) 


The  total  force  exerted  by  the  fluid  On  the  top 
plate  is  given  by  the  following  integration. 

fR 


F  =  j  (-%z)2.2H^T'^rdir 

•  Q 


(4.19) 


Substituting  equations  (4.10)/  (4.16)/  and  (4.13)  into 

(4.19)  and  performing  the  integration  one  obtains 

This  force  is  balanced  by  the  one  exerted  by  the  load/ 


r-  /  O 

F  =  m  (  3  -  J^) 


(4.21) 


Equating  equations  (4.20)  and  (4.21)  we  have 

rt 


^  ^ ^  (-  sr)  I J  ^  +  ^4  (-^J  “P  <e) 


(4.22) 
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Multiplying  this  equation  by  exp( -^ )  and  differentiating 
with  respect  to  time  gives  the  following  second  order 
dynamical  equation; 


1  ^  ^  r  I  J.  V  n  4.  ii  )1 


371VR: 


[l  +  f  4)"Jv  =  m3 


(4.23) 


When  R/H>>1,  equation  (4.23)  reduces  to 


fm  A,-^  +  [  1  f  ■^]  ^  +  f  ^ 


(4.24) 


where  /^(=  ^*^4  )  is  a  pseudo-fluid  density  corresponding  to 
the  amount  of  the  load  in  the  unlubricated  compressive  flow. 
The  group  arises  when  we  compare  the  load  inertia  to  the 
fluid  inertia,  that  is. 


load  inertia  =  m 


fluid  inertia  = 


H  dt 
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When  Ai  =  0,  equation  (4.24)  reduces  further  to 


P  A  4-  P  i  w  =  p  q 

dt*  '  dt  ^  2.  H*  ^  imC/ 


(4.25) 


This  aquation  is  to  be  compared  to  the  counterpart  for 
lubricated  flow,  equation  (3.29), 

,,  1  ^  o  iife  ^  X  /  _  4  m  3 

>°^i  ^  °  d-t  +-2^  4'  ■' 


TLR^ 


or  since  V=4H6b  from  equation  (3.2) 


Pi 


+  /‘-^  +  24^V  = 


(3.29A) 


I  ^  rn  H 

It  is  to  be  noted  that  is  equal  to  in  the 

lubricated  case. 


Equations  (4.25)  and  (3.29A)  have  very  similar  forms 
except  for  the  coefficients  of  the  third  terms.  Equation 
(4.25)  has  Vh*/  while  (3.29A)  has  ^/R^.  This  difference 
comes  from  the  different  origin  of  this  third  term  in  each 
equation;  the  term  in  equation  (4.25)  is  from  the  shear 
stress,  and  the  term  in  equation  (3.29A)  is  from  the  normal 
stresses . 

The  condition  under  which  equation  (4.24)  predicts 
the  occurrence  of  oscillatory  motion  is  given  by: 

Prn^  i  i  ^  )  “  4  /m  Ai  <  0 

If  the  fluid  inertia  is  included  in  the  more  general 
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development,  this  becomes 


/■*“(  i  + 


<  0 


r  ^  ^1  ^  /  I  -j-  A  JL^  ) 
b  >  (  I  +  i  /?»H*  / 


(4.26) 


/•  S  r  I  -t.  A  _?i  ^  ^ 


(4.27) 


Here  p*  =  /^  + 


Re  = 


HVP* 

1 


(modified  Reynolds  No.) 


Rd  = 


(Retardation  No.) 


De  = 


(Deborah  No.) 


When 


6 »(l  -f-  |-  -|^)^  r  equation  (4.24) 


following  approximate  oscillatory  solution, 


has  the 


The  oscillation  period  is  obtained  as  4/l.H 


or,  in 


P*Xi 

general,  as  4/lH^ 


This  again  has  a  different 


dependence  of  the  length  scale  from  the  period  in  the 
lubricated  case,  aR\  <■>,  • 


The  condition  (4.27)  is  conservative  for  the  reason 
that  has  been  discussed  in  Section  3.2.2;  that  is,  the  time 
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Numerical  solutions  of  equation  (4.24)  are  given  in 
this  section.  H/Hq  vs.  tiitie  curves  are  shown  in  Pig.  4.2  at 
various  values  of  relaxation  times  (  X|  =  0.003,  0.01,  0.03) 
and  zero  retardation  time,  and  compared  to  the  corresponding 
Newtonian  curve.  When  is  less  than  unity,  the  curve 

does  not  oscillate,  as  expected  from  the  condition  (4.29). 

As  0.6  increases,  the  oscillation  begins  and  becomes 

P*H 

severe.  The  oscillation  period  is  about  the  same  as  the 
calculated  value  of  T=47IHl/  .  Most  of  the  viscoelastic 
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curves  remain  below  the  corresponding  Newtonian  curve. 

The  effect  of  the  retardation  time  is  shown  in 
Fig.  4.3.  Hare,  the  relaxation  time  ( Ai )  is  kept  constant 
and  the  retardation  time  (Az)  varies  (  Az =  0.,  0-001,  0.003, 
0.01).  A  small  increase  of  the  retardation  time  reduces  the 
oscillation  amplitude  greatly,  as  expected  from  the 
condition  (4-30),  and  the  response  approaches  the  Newtonian 
curve  as  A^  approaches  Ai . 

4.2  Gontravariant  convected  Maxwell  fluid 
4.2.1  Governing  equations 
With  the  following  kinematics, 

Vr=  VrCk,  z,.t)  ^  \4  =  Vz(r,  z,t)  ,  =  0  (4.31) 

the  non-zero  components  of  the  deformation  rate  tensor  d 
are  6rr  >  <^96  •  ^zzf  <^rz>  from  which  we  also  assume  that 

the  non-zero  components  of  the  extra  stress  tensor  %  are 
'Trf  '799'  ^zz'  unlubricated  compressive  flow 

of  a  contravariant  convected  Maxwell  fluid  is  then  governed 
by  the  following  equations  and  the  boundary  conditions: 

Continuity ; 

T  +  Iz  =  0 


Momentum: 


(4.32) 


NEWTONIAN 


the  relaxation  time 


Constitutive : 


'?'e«  +A[^-27a,-^]=2^^ 


(4.35) 


(4.36) 


d3b  _2'7  -jy 

“of  9z  ^  arJ 

‘D7rz  ^  9Vz  y  aVr  - 
.  Dt  ar  az 

'fz\  ar 

(4.37) 


0  .... 

in  which  is  the  material  time  derivative 


Boundary  conditions: 


Vj.  =  0  at  z  =  0,  2H 


V,  =  0  at  z  =  0 


V,  =  -V  at  z  =  2H 


(4.39) 


=  0  at  r  =  0 


=0  at  r  =  0 


All  the  stresses  vanish  at 


the  free  surface 
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Here  we  consider  the  region  (0,R)x(0,2H)  (the  shaded  portion 
in  Fig.  4.1).  Note  that  we  do  not  have  flow  symmetry  with 
respect  to  z=H  (the  center  plane)  because  we  wish  to  include 
the  fluid  inertia  and  the  load  inertia,  while  the  top  plate 
is  moving  down  and  the  bottom  plate  is  stationary. 

In  the  linear  viscoelastic  case,  it  was  possible  to 
solve  for  the  stress  components  explicitly,  assuming  that 
the  velocity  profiles  are  the  same  as  for  a  Newtonian  fluid. 
Then,  substituting  into  the  momentum  equation,  we  were  able 
to  compute  the  closing  rate.  In  the  present  problem  we 
cannot  obtain  the  stress  components  explicitly  because  of 
the  nonlinearity  of  the  constitutive  equations.  Therefore, 
we  have  to  solve  seven  partial  differential  equations 
(4.32-4.33)  simultaneously  under  the  given  condition  (4.39). 
The  finite  element  numerical  technique,  which  has  been 
tested  in  Section  3.4,  will  be  applied  to  solve  this 
probl em . 

4.2.2  Finite  element  simulation 

In  section  3.4,  we  were;  able  to  solve  the  lubricated 
compressive  flow  of  a  contravariant  convected  Maxwell  fluid 
using  a  finite  element  numerical  technique.  We  now  wish  to 
solve  the  unlubricated  compressive  flow  problem  using  the 
same  finite  element  routine  but  with  different  boundary 
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conditions.  The  no-sIip  boundary  conditions  along  the  solid 
wall  (sea  (4.39))  are  imposed  in  the  present  problem,  while 
the  slip  boundary  conditions  were  used  in  the  lubricated 
case.  Since  the  no-slip  boundary  conditions  cause  the 
stress  singularity  at  the  edge  of  the  disk,  we  use  small 
elements  around  the  edge  of  the  disk  to  relax  the 
singularity  within  a  small  neighborhood  of  the  edge. 

The  initial  grid  used  in  the  calculation  and  the 
shapes  of  the  boundaries  at  later  times  are  shown  in 
Fig.  4.4  and  the  calculated  results  are  given  in  Fig.  4.5, 
in  which  H/Ho  is  plotted  against  time.  This  is-  compared  to 
the  linear  viscoelastic  case.  It  is  seen  from  this  figure 
that  the  overall  behavior  of  the  Maxwell  fluid  is  about  the 
same  as  in  the  linear  viscoelastic  case,  except  that  the 
Maxwell  case  shows  a  somewhat  smaller  amplitude  of  the 
oscillation,  which  may  be  due  to  the  non-linear  behavior  of 
the  Maxwell  fluid.  Both  viscoelastic  curves  remain  below 
the  corresponding  Newtonian  curve. 

The  shear  rate  in  this  flow  depends  upon  the 
position  and  the  time.  The  maximum  shear  rate  occurs  at  the 
edge  of  the  disk  and  at  the  time  when  the  closing  speed 
reaches  the  highest  value.  In  the  calculation  presented  in 
Fig.  4.4  and  4.5,  this  maximum  value  of  shear  rate,  Tln^x  »  is 
210  sec"'  at  time  t=0.22.  Since  the  relaxation  time( A ) 
used  in  this  calculation  is  0.1  sec,  A  turns  out  to  be 
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t  =  0.066 
H/Ho  =  0.866 


Fig.  4.4  The  initial  grid  and  the  shapes  of  the 
■boundaries  at  later  times  in  the  unlu'bricated 
compressive  flow  of  a  contravariant  convected 
Maxwell  fluid. 
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21.  This  high  value  of  ^max  A  ,  under  which  the  solution 
converges,  is  not  surprising  if  one  considers  that  this  flow 
is  an  oscillatory  transient  flow.  The  value  of  is 

zero  at  t=0,  it  increases  as  the  squeezing  goes  on,  and  it 
reaches  the  maximum,  then  it  decreases,  etc.  Higher  values 
of  were  tried,  but  the  solution  didn't  converge  after 

several  time  step  when  A  reaches  36.  In  the  steady 

state  calculation  of  Maxwell  fluid,  the  highest  value  of 
,  under  which  the  solution  converges,  has  been  known  as 
around  1  to  2  depending  upon  the  type  of  the  flow  ( see 
Crochet  and  Keunings,  1982;  Mendelson  et.  al . ,  1982; 

Viriyayuthakorn  and  Caswell, 1980) .  It  is  desirable  to  do 
the  computation  with  higher  j^oxA  in  the  future,  when  the 
convergence  problem  at  high  7\  t  which  is  one  of  the  major 
problem  in  the  numerical  calculation  of  viscoelastic 
materials,  is  resolved. 

At  this  point,  one  can  conclude  that  approximating 
the  overall  behavior  (H/Hq  vs.  time)  of  the  Maxwell  fluid  by 
the  linear  viscoelastic  prediction  is  favorable  at  the 
values  of  y^axA  at  least  up  to  21,  since  the  finite  element 
calculation  of  Maxwell  fluid  is  several  orders  of  magnitude 
more  expensive  than  th ?  linear  viscoelastic  case. 


CHAPTER  5 


EXPERIMENTS 

Experiments  on  the  unLubricated  compressive  flow  of 
Newtonian  and  viscoelastic  materials  under  a  constant  load 
have  been  carried  out  to  obtain  the  film  thickness  as  a 
function  of  time.  The  apparatus  and  the  materials  used  are 
described  in  Sections  5.1  and  5.2,  respectively.  The 
experimental  results  are  shown  and  compared  to  the 
theoretical  predictions  in  Section  5.3. 

Experiments  on  the  lubricated  compressive  flow  were 
attempted;  these  were  unsuccessful,  since  the  lubricant 
along  the  wall  surface  was  expelled  quickly  during  the  early 
moments  of  squeezing,  after  which  the  central  test  material 
began  to  stick  to  the  wall.  To  perform  this  type  of 
experiment  successfully,  the  viscosity  ratio  of  the 
lubricant  and  the  test  material  must  be  chosen  very 
carefully  (see  Section  3.5),  or  a  sophisticated  device  which 
can  supply  the  lubricant  to  maintain  almost  constant 
thickness  of  the  lubricant  layer  through  the  squeezings 
needs  to  be  designed. 

5 . 1  Apparatus 
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The  apparatus  used  in  the  experiments  is  shown  in 
the  picture  (Fig.  5.i).  There  are  three  parts;  the 
squeezing  equipment  {A),  the  thickness  measuring  device  (B), 
and  the  recording  equipment  (C) . 

5.1.1  Squeezing  equipment  (see  Fig.  5.2) 

The  squeezing  equipment  is  composed  of  a  rigid 
stationary  part  and  a  moving  part.  The  central  cylindrical 
rod  (Al),  at  the  bottom  of  which  a  flat  circular  disk  (A2) 
is  attached,  moves  through  two  linear  ball  bushings  (A3,A4) 
which  provide  the  straight  movement  of  the  moving  rod.  The 
test  material  is  placed  between  the  bottom  plate  (A5)  and 
the  circular  disk.  At  the  instant  t=0,  the  load  (all  the 
moving  parts)  on  the  material  is  released.  The  upper  plate 
falls  down  under  the  influence  of  the  gravitational 
acceleration,  with  the  test  material  being  squeezed  out. 

5.1.2  Thickness  measuring  device  (LVDT,  see 
Fig.  5.3) 

Continuous  measurement  of  the  film  thickness  during 
squeezing  is  accomplished  with  the  LVDT  (linear  variable 
differential  transducer),  which  is  manufactured  by  Schaevitz 
Eng.  Co.  The  specifications  of  the  LVDT  used  are  given  in 
Table  5.1.  The  LVDT  is  composed  of  two  parts,  the  body  ( B1 ) 
and  the  core  (S2).  The  body  is  fixed  on  the  stationary  part 
of  the  squeezing  equipment  and  the  core  is  attached  to  the 


Fig.  5.1  The  apparatus  used  in  the  squeezing 
experiments:  (A) squeezing  equipment,  (B) thickness 
measuring  device,  (C) recording  equipment. 


Fig.  5.2  The  squeezing  equipment:  (Al) cylindrical 
moving  rod,  (A2)Gircular  disk,  ( A3, A4) linear  ball 
bushings,  (a5) bottom  plate. 


Fig.  5.3  The  thickness  measuring  device 
(LVDT):  (Bl)body,  (B2)core. 
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Manufacturer  :  Schaevitz  Engineering  Co. 

Model  no.  ;  200  DC-D 

Calibration  data  by  manufacturer 

Linear  range  :  ±0.200  in. 

Sensitivity  :  51.500  V/in. 

Linearity  :  <0.25%  of  full  range 

AC  ripple  :  <10  mV  (max.) 

Dimensions 

Body  ;  0.75  in.  d  *  3.80  in.  1 

Core  :  0.187  in.  d  *  1.80  in.  1 

Weight 

Body  :  7  3  g 

Core  ;  5  g 

Operated  by  ±15  V  DC 


Table  5.1  The  specifications  of  LVDT  used 
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moving  cylindrical  rod,  so  that  the  core  moves  with  the  rod. 

5.1.3  Recording  equipment 

The  continuous  measurement  from  the  LVDT  is  sent  to 
the  oscilloscope  (Textronix  531A)  and  appears  as  a  moving 
spot  on  the  screen.  The  picture  of  the  trace  of  the  moving 
spot  is  taken  with  the  Tektronix  C-13  camera  and  a  Polaroid 
Land  pack  film  camera  back. 

5 . 2  Materials 

Two  Newtonian  fluids  have  been  used  as  the  standards 
to  -test  the  apparatus.  Three  different  viscoelastic 
materials  have  been  used  in  the  experiments  to  investigate 
viscoelastic  effects  in  the  compressive  flow. 

5.2.1  Newtonian  materials 

Viscasil  50000  :  This  is  a  viscous  silicone  fluid 
manufactured  by  General  Electric  Go.,  whose  viscosity  is 
constant  and  known  as  60000  cs  (at  25®C)  and  density  is  0.97 
g/cm^  . 

Dow  Corning  200  fluid,  12500  :  manufactured  by  Dow 
Corning  Co.  Its  viscosity  is  12500  cs  (at  25  G)  and  density 
is  0.975  g/cm^. 


5.2.2  Viscoelastic  materials 
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Silicone  Polymer  :  This  is  a  three  phase  material 
(silicone  resin,  plasticizer,  and  filler)  manufactured  by 
ICI,  England.  This  material  is  knov/n  to  have  rheological 
properties  close  to  those  of  linear  Maxwell  fluid,  based 
upon  the  oscillatory  shear  measurements  (see  also  Fig.  19-2 
of  Denn(1980) ) . 


The  oscillatory  shear  data  have  been  obtained 
through  the  courtesy  of  Dr.  K.  F.  Wissbrun  of  the  Celanese 
Research  Corporation.  Appendix  F.l  contains  a  tabulation  of 
the  storage  modulus(G'),  the  loss  modulus(G")  and  the 
absolute  value  of  the  complex  viscosity(  as  a  function  of 
the  circular  frequency (a>)  at  three  different  temperatures. 
By  taking  23°C  as  a  reference  temperature,  could  be 
superposed  on  to  one  mastercurve  with  the  help  of  the 
horizontal  shift  factors (^7)  (see  Gupta( 1980 ) ) ,  which  is 
shown  in  Fig.  5.4.  And  furthermore  we  know  from  the  linear 
viscoelastic  theory  that 


or 


(  I  + 


'^(jo  —  I  CO  (  I  +  ^ ^ 


(5.1) 


(5.2) 


Using  (5.2),  vs.  ci>dj  is  plotted  in  Fig.  5.5,  from  which 
one  can  see  that  the  viscosity  is  nearly  constant  in  the 
given  range  of  the  circular  frequency.  V7e  will  use  this 
constant  viscosity  to  analyze  the  squeeze  film  data  later. 
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TlA-227  :  A  concentrate  of  a  methacrylate  copolymer 
in  petroleum  oil,  manufactured  by  Texaco,  Inc.  The  steady 
shear  measurements  at  low  shear  rate  range  were  made  on  the 
Rheometrics  Mechanical  Spectrometer  (of  Celanese)  and  the 
capillary  measurements  were  made  at  high  shear  rates.  The 
shear  stress  and  the  normal  stress  data  are  tabulated  as  a 
function  of  the  shear  rate  in  Appendix  F.2.  The  end 
correction  (see  Bagley,  1957)  in  the  capillary  measurement 
was  unnecessary  since  the  capillary  tubes  used  were  long 
enough  (L/R=160,  268).  The  results  are  shown  in  Fig.  5.6  at 
various  temperatures.  At  27®C,  the  spectrometer  data  and 
the  capillary  data  agree  wall,  which  illustrates  the 
correctness  of  both  data.  The  material  shows  slightly  shear 
thinning  behavior ( n=0 . 86 )  at  high  shear  rates. 

PAA-water  solution  :  3.3  wt.  %  Separan  AP-30  (a 
partially  hydrolyzed  polyacrylamide  manufactured  by  Dow 
Chemical  Company)  in  water  solution  was  made  and 
characterized  with  a  Weissenberg  Rheogoniometer .  The  shear 
stress  and  the  first  normal  stress  data  are  given  in 
Appendix  F.3  and  shown  in  Pig.  5.7.  This  material  is  highly 
shear  thinning ( n=0. 267 )  and  highly  elastic. 

5 . 3  Experimental  results  and  discussion 

Experimental  results  on  Newtonian  fluids  are 
tabulated  in  Appendix  G.l  and  shown  in  Fig.  5.8,  in  which 
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Fig.  5*6  Shear  stress  vs.  shear  rate  of 
TLA-227 . 
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Fig.  5.7  The  shear  stress  and  the  first  normal 
stress  difference  of  3.3  wt.  %  PAA  in  water 
solution . 


Viscasil  60000 


122 


r  f  IS  f  iHo\  ts 

H/Ho  is  plotted  against  dimensionless  time,  t  y)  • 
They  agree  well  with  the  Stefan  eguation,  which  illustrates 
■the  quality  of  the  apparatus  and  technique  used. 


Experimental  results  on  viscoelastic  materials  are 
given  in  Appendix  G.2  and  shown  in  Figs.  5.9-5.13.  The 
corresponding  inelastic  curves  (analytical  or  numerical)  are 
also  shown  in  each  figure  for  comparison  purposes.  The 
values  of  for  each  run  are  also  given  in  Appendix  G.2. 

Figs.  5.9  and  5.10  show  the  squeezing  of  Silicone 
polymer.  It  is  seen  from  Fig.  5.9  that  the  material  seems 
to  be  squeezed  instantaneously  at  t=0.  The  details  of  this 
initial  movement  have  been  obtained  by  expanding  the  time 
axis  and  are  shown  in  Fig.  5.10,  in  which  the  linear 
viscoelastic  prediction  is  also  shown  for  comparison 
purposes.  Since  Tm«x  A  in  this  experiment  is  9,  we  expect 
that  the  Maxwell  prediction  would  be  close  to  the  linear 
viscoelastic  one  based  upon  the  conclusion  drawn  in  Section 
4.2.2.  It  shows  that  the  initial  movement  is  not  the 
instantaneous  squeezing,  but  the  oscillatory  squeezing,  and 
that  the  linear  viscoelastic  theory  predicts  the  correct 
oscillation  period,  but  a  larger  oscillation  amplitude  than 
the  experimental  result.  Both  experimental  results  on 
silicone  polymers  stay  below  the  corresponding  Newtonian 
curve  and  it  is  likely  that  this  material  is  close  to  a 
Maxwell  type  of  fluid,  as  expected  from  the  oscillatory 
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Fig.  5*9  H/HoVS.  time  in  the  unlubricated  squeezing 
of  silicone  polymer,  compared  to  the  corresponding 
Newtonian  fluid. 


124 


u 

d) 

(0 


•-I 

B 


Ld 


c 

•f-f  cn 

N  c 

d)  -H 

<D  4-> 
D  C!  !m 
0^0  0 
W  CU42 
CO  CO 
0) 

U  •• 

..  0 
u  0  CD 
CO 
d) 

x; 

4J 

O 
'H 
4J 
CO 

CD.  T3  ^ 
^  d)  ^ 
-P  CD 
ftl  O 

a  u 
S  CO 
o  -H 
CD  a  > 


»d 

d) 

4J 

rd 


u 

f— ( 

d 

d 


CO 


rd 

o 


o 

-p 


c 

•H 


s 

•H 

+J 


-  u 

u  id 

d)  CD 


0 

C 

C! 

\  d)  Id 

c 

0 


d) 

CO 

d 

0 

cu 


o  O  d3  CO 


rH  -H 


lO  *H 
•  <0 
CD 


•H 

d 

o 

aJ 

y 


CD 

^  B 

-H  CD  -H 

0  S  4-> 


OH/H 


126 


OH/H 


of  TLA.-2  27,  compared  to  the  corresponding  power-law  case. 
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shear  data. 


The  squeezing  curves  of  TL?i-227  are  shown  in 

Figs.  5.11  through  5.13.  One  can  see  that  the  experimental 
curves  stay  above  the  corresponding  inelastic  curves. 
Fig.  5.14  shows  the  result  for  the  PA\-water  solution,  which 
has  the  similar  characteristic  as  TL\-227.  Most  data  on 
these  two  solutions  show  an  inflection  point  which  probably 
reflects  a  very  weak  oscillation.  In  Section  4.2.2  we  have 
seen  that  the  Maxwell  model  predicts  the  faster  squeezing  at 
(=21)  which  is  higher  than  those(=2 10)  in  these 
experiments.  Thus,  it  is  not  likely,  at  least  in  this  type 
of  transient  flow,  that  these  two  polymer  solutions  are 
close  to  Maxwell,  or  in  general  White-Metzner ,  type  of 
fluids . 


It  is  desirable  to  try  more  realistic  and  thus 
complex  models,  which  can  describe  the  transient  behavior 
more  precisely,  in  finite  element  calculations  in  the 
future.  The  Phan-Thien-Tanner ,  Johnson-Segalman,  and 
Marrucci  structural  models  would  all  be  candidates  for  this 


purpose . 


CHAPTSR  6 


SUMMARY  AND  RECOMMENDATIONS 

6 . 1  Summary  of  present  work 

1.  The  unlubricated  compressive  flows  of  Newtonian 
fluids  and  power-law  fluids  have  been  simulated  by  a  finite 
element  technique,  verifying  that  the  Stefan  equation 
(Newtonian)  and  the  Scott  equation  (power-law)  based  upon  an 
assumption  of  parallel  squeezing  are  good  provided  the  R/H 
ratio  is  large  enough.  When  R/H  is  small,  the  edge  effect 
causes  slower  squeezing  than  predicted  from  the  Stefan 
equation  or  the  Scott  equation. 

2.  Xfhen  there  exists  a  substantial  transverse 
viscosity  gradient  in  the  fluid  charge,  two  different  flow 
regimes  are  predicted  depending  upon  whether  the 
dimensionless  group,  3  (=  vr^-^)  is  small  or  large  compared 
to  unity.  When  S  is  small  compared  to  unity,  the  parallel 
squeezing  assumption  is  valid  and  the  maximum  velocity 
occurs  at  the  center  plane.  When  S  is  large  compared  to 
unity,  the  parallel  squeezing  assumption  breaks  down  and  the 
maximum  velocity  occurs  in  the  low  viscosity  fluid  near  the 
disks.  Thus,  a  new  analysis  is  necessary  in  this  case.  In 
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the  two  fluids  case,  new  assumptions  other  than  the  parallel 
squeezing  assumption  have  been  made  to  derive  an  analytical 
solution  which  is  found  to  be  in  good  agreement  with  the 
finite  element  numerical  solution. 

3.  The  partially-filled  compressive  flow  of 
Newtonian  fluid  has  been  solved  -numerically.  The  flow 
patterns  in  the  partially-filled  compressive  flow,  or  the 
flow  in  the  mold  cavity,  are  essentially  the  same  as  those 
one  observes  in  the  fully-filled  case  except  near  the  front 
:  here  we  observe  the  expected  fountain  flow  phenomenon  in 
the  partially-filled  case. 

4.  The  lubricated  compressive  flow  of  linear 
viscoelastic  material  has  been  solved  analytically  and 
numerically.  The  compressive  motion  under  a  constant  amount 
of  load  may  or  may  not  be  oscillatory,  depending  upon  the 
conditions  as  predicted  by  equations  (3.49-52).  Ifhen 
oscillation  occurs,  it  is  due  to  the  combined  effects  of 
inertia,  the  elasticity,  and  the  viscosity.  The  oscillation 

,rpT 

period  is  given  by  TtRy  .  The  retardation  time  plays  an 
important  role  in  the  damping  of  the  oscillatory  motion. 

5.  The  lubricated  compressive  flows  of  non-linear 
viscoelastic  model  fluids,  including  the  contravariant 
convected  Maxwell  model,  Vfliite-Metzner  model, 
Johnson-Segalman  model  (with  and  without  retardation  time) , 
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and  Marrucci  structural  model,  have  been  solved  numerically. 
The  contravariant  convected  Maxwell  model  and  White-Metzner 
model  predict  faster  oscillatory  squeezing  than  the 
corresponding  inelastic  cases,  except  at  very  small  H/Hq 
under  the  extraordinarily  high  loading  conditions.  The 
Johnson-Segalman  model  and  the  Marrucci  structural  model 
predict  slower  oscillatory  squeezing  than  the  corresponding 
inelastic  cases  at  high  squeezing  speed.  These  different 
responses  seem  to  be  due  to  the  stress  overshoot  phenomena, 
since  Johnson-Segalman  model  and  Marrucci  structural  model 
do  predict  stress  overshoot,  but  the  contravariant  convected 
Maxwell  model  and  V'Hiite-Metzner  model  do  not. 


6.  The  unlubricated  compressive  flow  of  linear 
viscoelastic  materials  has  been  solved  analytically, 


assuming 

that 

the  velocity  field 

is 

the  same 

as 

the 

Newtonian 

one . 

The  behavior  is 

very 

similar 

to 

fhe 

lubricated 

case 

except  that  the  principal 

length 

scale 

is 

different  from  the  lubricated  case.  When  oscillation 

JFx 

occurs,  the  oscillation  period  is  given  by  6"^  • 


7.  The  unlubricated  compressive  flow  of  the 
contravariant  convected  Maxwell  model  has  been  solved 
numerically,  using  a  finite  element  technique.  The  results 
are  about  the  same  as  in  the  linear  viscoelastic  case  except 
that  the  oscillation  amplitude  of  the  Maxwell  fluid  is 
somewhat  smaller  than  for  the  linear  viscoelastic  material. 
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■  i^n  -hTie  visGOQl^stic 

3.  Squeezing  experiments  on  tne 

.atarials  have  ehcwa  two  aifferent  Oehaviots,  silicone 

polymer,  which  can  be  well  characterlrea  by  the  linear 

Maxwell  model,  shows  faster  oscillatory  squeezing  than  the 

corresponding  Newtonian  cases,  while  TL!l-227  and  flVU-water 

solution  Show  weaK  bounces  and  slower  squeesing  than  the 

corresponding  inelastic  cases.  The  behavior  of  silicone 

polymer  is  predictable  from  the  Maxwell  model,  but  the 

responses  of  TL.4-227  and  P!iA-watar  solution  are  different 

from  those  predicted  by  the  Maxwell  model.  To  explain  the 

latter  behavior,  we  might  need  a  constitutive  equation  which 
can  predict  the  transient  responses,  such  as  the  stress 

overshoot,  xnora  precisely* 

9,  in  conclusion,  the  slower  squeeslng  of 
viscoelastic  materials,  or  higher  load-bearing  capacity  of 
viscoelastic  lubricants,  seams  to  be  not  due  to  the 
elasticity  of  the  material,  but  due  to  the  transient 


responses  of  the  material 


In  other  words,  the 


dimensionless  group  alone  may  not  be  the  proper 

dimensionless  group  to  describe  the  whole  features  of  the 
squeezing  of  viscoelastic  materials.  We  require  the  ratio 
of  Deborah  number  to  Reynolds  number  to  describe  the 
csclllatory  behavior  and  also  require  at  least  one  more 
dimensionless  group  which  can  describe  the  transL.-nt 
responses  of  the  materials  such  as  stress  overshoot 


phenomenon . 
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6 . 2  Recommendations  for  future  work 

1.  The  use  of  the  Phan-Thien-Tanner  model,  Johnson- 
Segalman  model,  or  Marrucci  structural  modal  is  recommended 
for  the  further  study  of  the  compressive  flow  problems. 

2.  In  all  transient  problems  it  may  be  important  to 
use  constitutive  aquations  which  describe  the  transient 
fluid  behavior  precisely. 

3.  The  development  of  a  constitutive  equation,  which 
can  describe  the  transient  responses  more  precisely,  is 


desirable . 
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APPENDIX  B 


GOVERNING  EQUATIONS  IN  THE  LUBRICATED  GPMPRESSIVE 
FLOv^  OF  NON-LINEAR  MODEL  FLUIDS. 

Assuming  the  kinematics  given  by  equations  (3.1-3), 
the  lubricated  compressive  flow  of  non-linear  model  fluids 
are  governed  by  equations  (3.19,  16,  22),  and  the 

corresponding  constitutive  equations  which  calculate  Ter  and 
Tzz .  The  constitutive  equations  which  will  be  considered 
here  include  White— Metzner  model,  Johnson— Segalman  model 
(with  or  without  retardation  term),  and  structural  model. 
Since  the  contravariant  cbnvected  Maxwell  model  is  just  a 
special  case  of  White-Metzner  model  (when  (]l^)='^o)' 
Johnson-Segalman  model  (when  3=0)/  it  will  not  be  dealt  with 
separately. 


B.l  White-Metzner  model 


Rewriting  equations  (3.19,  16,  22), 

d&h  W  ( 5 1  6 H )  ~  -^TIP t  ^ R  H  ^  -  TlR  (T!r~ ^z) 

dt  = 


dH 

dt 


=  -2H% 


(3.19) 


(3.16) 


d  —  D  g, 

TT  - 


(3.22) 


The  Trr  and  7zz  in  (3.19)  are  calculated  using  White-Metzner 
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constitutive  equation  given  by  equation  (3.6  ) 

Z  +  ^  2  =  <3-6  > 

In  the  lubricated  compressive  flow, 

/  6j,  0  0  ^ 

d  =  VV  +  VV"^ )  =  0  4  0  ( B .  1 ) 

^0  0  -2.4  y 


and 

Ij  =-'j[(trdf  -  tr(d^)]  =  +3£i,  (B.2) 

If  one  assumes  that  the  viscosity  is  given  by  the  power-law 
relationship,  that  is, 

n-i 

(Id)  =  k|+4]Ij|  ^  {B.3) 

where  K  is  the  consistency  factor  and  n  is  the  power-law 
index,  the  radial  and  axial  components  of  equation  (3.65) 
are  written  as 


2<5-4- 


+  ^4  7rr 


d  '^2Z 
dt 


~  4  '^ZL 


(B.4) 

(B.5) 


Equations  (3.19,  16,  22)  and  (B.4,  5)  are  all  first 
order  ordinary  differential  equations  and  they  can  b<=‘  solved 
simultaneously  using  the  Runge-Kutta  method  or  Gear's 


method . 


B. 2  Johnson-Sagalman  model 

Equations  (3.19,  16,  22)  remain  the  same  and  the 
stress  components  (rr  and  zz)  of  Johnson-Segalman 
constitutive  equation  with  retardation  term  (3.55),  in  the 
lubricated  compressive  flow,  are  given  by 

7rr<-  a,[4^  -20-J)4r,r]  = 

7=  +  A,  L^+4(i-w4  T=]  =-«J[4+A4^  +  'K/-J)6/}] 

or 

Again  it  is  not  difficult  to  solve  the  equations  (3.19,  16, 

22)  and  (B.6,  7)  simul tanuously . 

B.3  Structural  model 


First  we  need  to  rearrange  equation  (3.19)  in  the 
following  form. 


since  the  constitutive  equations  are  written  in  terms  of 
Tj-(.-(=  'Ti-r/G)  and  '<zz(=  'Zcz/G)  . 

In  the  lubricated  compressive  flow,  the  constitutive 


equations  (3.67)  become 
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'^rr  +  -  ^ehin-l  = 

7z^  fXL^  +  4€i,  7^j  = 

^  X  ^  A  =  -Ae>  -X 

dx  '  /,  1  ^  '  /IT  (B.IO) 

^=  — (/-x) -ax-^ 

All  the  equations  are  the  first  order  ordinary  differential 
equations  again,  which  can  be  solved  easily. 

Computer  programs  for  each  model  are  found  in  the 


Appendix  E.3.  Here,  Gear's  method  has  been  used  to  solve 
given  ordinary  differential  equations  simultaneously. 


APPENDIX  C 


THE  LUBRICATED  COMPRESSIVE  PLOW  OP  MAXWELL 
PLUID  UNDER  VERY  HIGH  LOADING  CONDITIONS 

The  numerical  results  of  the  lubricated  compressive 
flow  of  Maxwell  fluid  under  very  high  loading  conditions  are 
given  here.  The  geometry  and  the  material  properties  are 
kept  constant  as  given  below.  The  load  varies  from  100  to 
10000. 

Geometry  and  material  properties ' : 

R  =  3 . 

Ho  =  0.05 

/  =  1. 

^  =  10. 

A  =  0.1 
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1.  m=100. 


(H/Ho)Max 


t 

(H/Hq  ) 

(H/Ho)]yj^xwall 

» 

0.00 

1.000 

1.000 

0.00 

1.00 

0.01 

0.899 

0.893 

10.93 

0.99 

0.02 

0.694 

0.665 

17.55 

0.96 

0.03 

0.507 

0.459 

18.89 

0.91 

0.04 

0.327 

0.318 

17.59 

0.85 

0.05 

0.281 

0.228 

15.64 

0.81 

0.07 

0.175 

0.131 

12.13 

0.75 

0.1 

0.109 

0.710E-01 

8.60 

0. 65 

0.15 

0.522E-01 

0.358E-01 

5.46 

0.69 

0.2 

0.323E-01 

0.226E-01 

3.86 

0.70 

0.3 

0.161E-01 

0.124E-01 

2.36 

0.77 

0.4 

0.979E-02 

0.830E-02 

1.72 

0.85 

0.6 

0.473E-02 

0.474E-02 

1.18 

0.99 

0.8 

0.286E-02 

0.313E-02 

0.93 

1.09 

1.0 

0.191E-02 

0. 230E-02 

0.79 

1.20 

m=500 . 

• 

(H/Ho)Max 

t 

Newt 

(H/Ho)Max 

^^b^Max 

( H/Hq ) Newt 

0.00 

1.000 

1.000 

0.00 

1.00 

0.01 

0.746 

0.  738 

30.13 

0.99 

0.02 

0.368 

0.345 

40.18 

0.94 

Q.03 

0.187 

0.167 

32.12 

0.89 

0.05 

0. 721E-01 

0.604E-01 

20.09 

0.84 

0.1 

0.188E-01 

0.155E-01 

9.48 

0.82 

0.2 

0.488E-02 

0.449E-02 

4.18 

0.92 

0.3 

0.221E-02 

0.233E-02 

2.61 

1.05 

0.4 

0.126E-02 

0.149E-02 

1.93 

1.18 

0. 6 

0.569E-03 

0.793E-03 

1.33 

1.53 

0.8 

0.324E-03 

0.497E-03 

1.04 

1.53 

1.0 

0.209E-03 

0.341E-03 

0.86 

1.63 
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m=iOOO . 

(H/Ho)Max 

t 

^  ^  Max 

0.00 

1.000 

1.000 

0.00 

1.00 

0.01 

0.677 

0.670 

41.76 

0.99 

0.02 

0.251 

0.236 

51.40 

0.94 

0.03 

0.110 

0.988E-01 

36.51 

0.90 

0.05 

0.376E-01 

0.327E-01 

21.26 

0.87 

0.1 

0.903E-02 

0.796E-02 

9.72 

0.88 

0.15 

0.398E-02 

0.372E-02 

6.02 

0.93 

0.2 

0.223E-02 

0.225E-02 

4.27 

1.01 

0.4 

0.560E-03 

0.728E-03 

1.98 

1.30 

0.6 

0.251E-03 

0.381E-03 

1.36 

1.52 

0.8 

0.142E-03 

0.236E-03 

1.06 

1 . 66 

1.0 

0.911E-04 

0.161E-03 

0.87 

1.77 

m=10000 . 

(H/Ho)Max 

t 

(H/Ho)uewt 

(H/HQ)j^ax 

( ^  Max 

0.00 

1.000 

1.000 

0.00 

1.00 

0.01 

0.541 

0.539 

80.30 

0.996 

0.02 

0.552E-01 

0.534E-01 

87.82 

0.97 

0.03 

0.156E-01 

0.150E-01 

47.25 

0.96 

0.05 

0.414E-02 

0.399E-02 

24.00 

0.96 

0.07 

0.187E-02 

0.184E-02 

15.85 

0.98 

0.1 

0.846E-03 

0.859E-03 

10.30 

1.02 

0.15 

0.353E-03 

0.386E-03 

6.31 

1.09 

0.2 

0.193E-03 

0.228E-03 

4.46 

1.18 

0.4 

0.465E-04 

0.703E-04 

2.06 

1.51 

0.6 

0.204E-04 

0.358E-04 

1.41 

1.75 

0.8 

0.114E-04 

0.219E-04 

1.09 

1.92 

1.0 

0.729E-05 

0.147E-04 

0.89 

2.02 

APPENDIX  D 


FINITE  ELEMENT  FORMULATIONS. 

Since  our  flow  of  interest  is  the  ax i symmetric  flow 
(in  r-z  coordinate),  we  will  limit  ourselves  to  the 
axisymmetric  flow  of  incompressible  fluids.  First,  we  will 
consider  the  flow  of  generalized  Newtonian  fluid.  The  flow 
of  viscoelastic  materials  will  then  be  considered.  The 
excellent  treatments  on  these  subject  are  given  by 
Crochet( 1981 ) ,  and  Zienkiewicz{ 1977 ) . 

D.l  Generalized  Newtonian  flow 


Governing  equations  :  In  the  compressive  flow  of 
inelastic  generalized  Newtonian  fluid,  the  fluid  inertia  is 
considered  negligible  (Re<<l).  The  momentum  equations  then 
become  the  quasi-steady  state  equation,  even  though  the  flow 
itself  is  time  dependent.  The  momentiAm  equations  are  given 
by 


iL 

dr 


•h 


+ 


T  ~ 


r  9z 


-i-  ff,  =  0 


d7^ 

dZ 


•h  -f,  -  0 


(D.l) 


(D.2) 


where  f^  and  f^  represent  the  components  of  the  body  force 
per  unit  volume. 
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The  constitutive  relations  of  generalized  Newtonian 


fluids  are  given  by 


~  2  "^  ( ]ij  )  d,,^ 


(D.3) 


=  2  I  {  id  )  dj 


(D.4) 


72Z.  =  2  (Id  )  dzz 


(D.5) 


Trz  ~  2  ( Id  ^ 


{D.6) 


wnere 


(tr  d)  -  tr( 


5^) 


(  d^  +  dga  +  dzz  )  +  d(.z 


(D.7) 


^lrr~ 


(D.8) 


dao  — 


(D.9) 


dzz  - 


(D.IO) 


^rz-  Z  [  3Z  ^  dr 


(D.ll) 


The  mass  conservation  is  given  by 


I  ^  ^  ?.Vk.  =  0 
9r  ^  r  ^  az  ^ 


(D.12) 


Galerkin  formulation  :  Let  us  consider  the 


approximations  of  ,  and  p  given  by 
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I 


“it- 

p  =  i 


(D.13) 


(D.14) 


(D.15) 


where  'ik-  and  are  the  interpolation  functions  for  the 
velocity  components  and  the  pressure,  respectively,  and  t^-, 
Vj  /  and  Pj  are  the  nodal  values  of  the  velocity  components 
and  the  pressure.  The  Galerkin  formulation  of  the  momentum 


equations  is  then  given  by 


<<-f;  0  (0.16) 


<  ,  ~  sS:  T  + 


+  fz  y  -  0 


(D.17) 


where  the  brackets  denote  the  integrations  over  the  given 
domain . 


Performing  the  integration  by  parts  on  equations 
(D.16)  and  (D.17)  and  using  the  divergence  theorem,  one 


obtains 


<  '■fr,  +  +  <  I-# ,  -Jk  >  +  <ti  ,%,>-<  ll-i ,  P> 

=  <rfi,  fr>  +  «rl!'i  , 


where  t^.  and  tj  are  the  r  and  z  component  of  the  contact 
force  vector  and  the  double  brackets  denote  the  integrations 
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along  the  boundary  of  the  domain. 


Substituting  the  cohstitutive  relations  (D.3-6)  into 
equations  (D.18)  and  (D.19)  gives  the  final  form  of  Galerkin 
equations, 


-<^i,  P>  =  <  ri^i^  f^>  +  « 

=  <  ■Pz>  +  <  ±2  » 


(D.20) 


(D.21) 


Galerkin  form  of  the  mass  conservation  becomes 

✓  Jv  4_  4.  V  _  o 


y  1  e/vr  ,  Vr  .  <iV3 

<r-f‘,  ^  TT  +  T  az 


(D.22) 


Further  substituting  the  approximations  of  (D. 13-15) 
into  (D. 20-22),  one  obtains  the  following  algebraic  system 

I  (  AijUj  +  -  i  D^-  Pj  =  X;  ^  i:|~M  (D.23) 


-11 

J  J-. 


O'** 


i  :  1  ~  M 

(D.23 ) 

i  :  i 

(D.24) 

/  :  1 

(D.25) 

where 


Aij  ~  <  7  TF,  Tr  >  ^  ^  az  ,  az  ^  ^  ^  7  ;  r  ^ 

R  -  2.<^r^  lty+^o,rdL  i^> 

^b-<^4,  h> 

Xi  ~  <  m/>j  ^  >  +  «  n/';,  ^  Yi  =  <  >"V';  j  "^z  » 


152 


The  solution  of  the  system  (D. 23-25)  can  be  obtained 
by  a  certain  iteration  technique,  since  the  viscosity  is  a 
function  of  the  velocity  field  in  general.  One  can  assume  a 
constant  viscosity  to  obtain  the  velocity  field  in  the  first 
iteration  and  the  new  viscosities  are  obtained  from  the 
previously  obtained  velocity  field.  One  then  repeats  this 
procedure  until  the  solution  or  the  viscosity  converges 
within  a  given  error  allowance.  This  technique  has  been 
used  successfully  to  solve  the  unlubricated  compressive 
flows  of  the  power-law  fluids  in  Section  2.3. 

solving  equations  (D. 23-25)  requires  the  calculation 
of  all  the  matrices,  which  are  the  integrations  in  the 
domain,  and  the  selection  of  the  appropriate  elements  and 
shape  functions  i'l  and  .  Numerical  integration,  using 
quadrature  points,  is  the  effective  way  of  performing 
integations.  Two  types  of  elements  have  been  successful  in 
the  past.  The  first  element  is  a  triangle  on  which  the 
velocity  components  are  represented  by  complete  second  order 
polynomials,  while  the  pressure  is  given  by  complete  first 
order  polynomials.  The  second  element  is  a  quadrilateral 
with  biquadratic  velocity  components  and  bilinear  pressure. 
The  functions  and  4*',  are  given  in  Tables  D.l  and  D.2  for 
the  triangular  and  the  quadrilateral  element,  respectively. 
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D.2  Viscoelastic  flow  ;  Maxwell  fluid 

Governing  equations  :  Let  us  consider  the 
compressive  flow  of  the  contravariant  convected  Maxwell 
fluid.  One  then  requires  to  retain  the  material  time 
derivative  terms  in  the  momentum  equations  and  in  the 
constitutive  equations.  They  are  given  in  the  following. 


Momentum ; 


^  D-t 


££.  4.  _L  J-Jry  )  -  4-  f 

.  I  d  I^'T  \  ,  .  r 


(D.26) 


(D.27) 


Constitutive : 


PT^ 

Dt 

Z7  ^ 

^'7  ^'^'•1 

az  J 

=  3. 

/  ar 

(D. 

.28) 

Dt 

_9  7  1 

r 

(D. 

.29) 

A[ 

PTa 

Ot 

—2T 

-i'izz 

aVz  T 
az  J 

=  2  ■ 

/ 

(D. 

.30) 

■A[ 

D7w 

Dt 

y 

■'“■az 

.-7„ 

\ 

(D. 

1 

.31) 

+  ^1 

The  continuity  equation  is  given  by 


aK  ^  r  ^  az  ^ 


(D.32) 


A  mixed  finite  element  method  ;  In  this  method, 
first  proposed  by  Kawahara  and  Takeuchi,  the  extra  stress 
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components  are  considered  as  unknown  functions  to  the  same 
degree  as  the  velocity  components  and  the  pressure.  Thus, 
we  consider  7  unknown  fields  for  axisymmetric  flow;  4 
stress  components,  2  velocity  components,  and  the  pressure. 
We  use  the  following  approximations  for  these  unknown 


fields , 

(D.33) 

Vz  ~  (D.34) 

p  X  Pj  (D.35) 

(D.35) 

7zz-2.Sj7j  (D.37) 

(D.38) 

(D.39) 


The  Galerkin  form  of  the  constitutive  and  field  equations 
may  then  be  obtained  as  follows. 


<  rr,, 


-27 

2- (ft 


'k] 


(D.40) 


<  r7;  ^  Ofc*  r" ] 


(D,41) 


<  PTi 


2  =  0 


(D.42) 
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+  >  =  0 

(D.44) 

,  P>  =  ^«rii]±r» 


I  \  cfr  / 

7t7:>  +<t,  7,6 

-<i'\,  P>  =  <'-Vj,  f^>  +«1-V';^  tf 


<U-  ^ 


#>+<'-1^,  -^> 


1.45) 


By  replacing  V^  ,  y*.  /  '7’rr  /  'Tiz/  7fz.  7^9  .  and  p  in  terms  of 

their  approximations  given  by  (D. 33-39)  it  is  possible  to 
obtain  a  non-linear  algebraic  system  of  equations  in  terms 
of  the  unknowns,  Uj  ,  Vj  ,  Rj  ,  Sj  ,  Tj  ,  Qj  ,  and  ,  which  is 
solved  by  means  of  Newton-Raphson  iterative  method. 


The  material  time  derivatives  are  treated  by  a 
finite  difference  scheme  in  the  time  coordinate.  Their 
discretization,  using  the  implicit  three  point  recurrence 
scheme  with  variable  time  steps,  is  given  by  equation 
(3.72). 


The  same  types  of  elements  as  used  in  the 
generalized  Newtonian  flow  have  been  shown  to  work  well  with 
the  mixed  method,  using  the  same  interpolation  functions  for 
the  velccity  components  and  the  stress  components;  that  is, 

7i  ■ 


parent  triangular  element 

s 

4  =  ^ 

i/,  =  I  -3(5+3)  tiff+^j"' 

A  =  J(-25-l ) 

%  = 

A  =  4S(i-3-J) 


Table  D.l  Shape  functions,  and 
the  parent  triangular  element. 


parent  quadrilateral  element 


=  5n+l)^n+J)/4  ,  ii  =  ■'S(l-S)y\  +  J)/4 

^3=  5f/-5)^n-^)/4  3  A  =  -i(/+5;7(:i-)V4 

Table  D.2  Shape  functions,  and  ,  in 
the  quadrilateral  element. 


APPENDIX  S 


COMPUTER  PROGRAMS 

The  computer  programs  for  the  lubricated  compressive 
flows  of  the  White-Metzner,  Johnson-Segalman ,  Marrucci 
structural  models,  and  their  corresponding  inelastic  cases 
are  given  here.  Subroutine  DGEAR,  the  differential  equation 
solver  by  Gear  method,  from  IMSL  subroutine  package  has  been 
used  in  these  programs.  Each  program  deals  with  the 
following  case; 

Program  Description 

WM3LP.F0R  White-Metzner  model  with 

power-law  viscosity 

PWL3LP.P0R  inelastic  power-law  fluid 

JS3LP.P0R  Johnson-Segalman  model 

NJSSLP.FOR  inelastic  case  corresponding  to 

Johnson-Segalman  model 

SM3LP.F0R  Marrucci  structural  model 

NSMSLP.FOR  inelastic  case  corresponding  to 


Marrucci  model 
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Programs  used  in  the  lubricated  compressive  flows 


C**** ***************************************** ***************** 

C 

C  PROGRAM  WMSLP.FOR 

C 

c  LUBRICATED  SQUEEZING  OF  WHITE-METZNER  FLUIDS. 

C 

C*********** ************************************************** 

INTEGER  N,METH,MITER, INDEX, IWK{5) ,IER,IRV 

REAL  Y(5) ,WK(115) ,T,TOL,TEND,DELT 

DOUBLE  PRECISION  PRINT, THHO 

EXTERNAL  FCN,FCNJ 

COMMON/AB/A,B 

COMMON/AM/AMASS 

COMMON/G/GM 

COMMON/VR/VI SC , GSHEAR , PLI 

Q* *********************** ************* 

TYPE  10 

10  FORMAT ('  INITIAL  CONDITIONS'/’  ENTER  T0,EB0,PRR0 
1,PZZO,HO,RO' ) 

ACCEPT  *,T0,EB0,PRR0,PZZ0,H0,R0 
TYPE  11 

11  FORMATC  ENTER  DENSITY,  SHEAR  MODULUS,  VISCOSITY  OR' 

1  '  CONSISTENCY  FACTOR,  POWER  LAW  INDEX,  MASS') 

ACCEPT  * , DENS , GSHEAR , CONF , PLI , AMASS 
VISC=CONF*(12. )**( (PLI-1. )/2. ) 

TYPE  12 

12  FORMATC  ENTER  TOLERANCE,  DELT' ) 

ACCEPT  *,TOL,DELT 

TYPE  15 

15  FORMATC  ENTER  NPRT, INTERVAL' ) 

ACCEPT  *,NPRT,AINTR 

TYPE  16 

16  FORMATC  ENTER  FILENAMES  FOR  PRINT  /  T  VS.  HHO' ) 

ACCEPT  117, PRINT, THHO 

117  FORMAT (AlO/AlO) 

Q* *********************************** 

OPEN ( UNI T= 31, DEVI CE='DSK' ,FILE=PRINT) 

OPEN ( UNI T*3  2 , DEVI CE» ' DSK ' , FI LE»THH0 ) 

WRITE (31, 17) 

17  FORMATC  LUBRICATED  SQUEEZING  OF  WHITE-METZNER  FLUID') 
WRITE (31, 25)  DENS, GSHEAR, CONF, PLI , AMASS 

25  FORMATC  DENSITY  =  ’,E14.5/'  SHEAR  MODULUS  =  ’,E14.5/ 

1  '  CONSISTENCY  FACTOR  =  ',E14,5/'  POWER  LAW  INDEX  =  ', 

1  E14.5/'  MASS  =  ' ,E14.5) 

WRITE (31, 27)  TOL 

27  FORMATC  TOLERANCE  =  ',E14.5) 

^** *********************************** ****** 

PI=3. 14159265357989 
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N=5 

METH*1 

MITER=0-- 

INDEX=1 

WRITE (31, 28)  METH, MI TER, INDEX 

28  FORMAT(/'  METH  =  ',14/’  MITER  =  ',14/'  INDEX  =  ',14/) 


T=T0 

Y{1)=EB0 

Y(2)=PRR0 

Y(3)=PZZ0 

Y(4)=H0 

Y(5)=R0 

GM= 980.* AMASS 

A=PI*DENS 

B=PI*VISC 

EB=Y(1) 

PRR=Y(2) 

PZZ=Y(3) 

H^Y(4) 

R=Y(5) 

HH0=H/H0 
WRITE (5, 20) 

WRITE (31, 20) 

20  FORMAT(//7X, 'TIME' ,12X, 'EB' ,11X, 'PRR' ,11X, 'PZZ' , 

1  13X,'H' ,13X,'R' ,10X,'H/H0' , lOX, ' DELT’ /) 

WRITE(5,21)  T, EB, PRR, PZZ, H,R,HHO, DELT 
WRITE(31,21)  T, EB, PRR, PZZ, H,R,HHO, DELT 
WRITE(32,*)  T,HHO 

21  FORMAT(8E14.5) 

C********** ************************** 

DO  100  I=1,NPRT 
TEND=FLOAT( I ) *AINTR+T0 

CALL  DGEAR ( N, FCN, FCNJ,T, DELT, Y, TEND, TOL, METH, MI TER, 
1  INDEX, IWK,WK,IER) 

EB=Y(1) 

PRR=Y ( 2 ) 

PZZ*Y(3) 

H=Y(4) 

R=Y(5) 

HH0=H/H0 

WRITE (31, 21)  TEND, EB, PRR, PZZ, H,R,HH0, DELT 
WRITE(32,*)  TEND,HH0 

100  WRITE(5,21)  TEND, EB, PRR, PZZ, H,R,HH0, DELT 

CLOSE(UNIT=31,DISPOSE='PRINT' ) 

CLOSE ( UNI T=32) 

STOP 

END 

Q*  ********************************************  * 
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SUBROUTINE  FCN(N,T, Y, YPRIME) 

Q********************************************* 

INTEGER  N 

REAL  Y{5),YPRIME(5) ,T 
COMMON/AB/AfB 
COMMON/AM/AMASS 
GOMMON/G/GM 

COMMON/VR/VI SC , GSHEAR , PLI 

Q* ************************** *********** 

PI=3. 14159265357989 
EB=Y{1) 

EB2=EB*EB 

PRR=Y(2) 

PZZ=Y(3) 

H=Y(4) 

H2=H*H 

R=Y(5) 

R2=R*R 

R4=R2*R2 

YPRIME ( 1 ) = (GM+4 . *AMASS*H*EB2-A*R4*EB2/4 . 

1  +4.*A*R2*H2*EB2/3.-PI*R2*(PRR-PZZ) ) 

1  /( A*R4/4 .+2 . *A*R2*H2/3 . +2 . *AMASS*H) 

YPRIME ( 2 ) =2 . *GSHEAR*EB-GSHEAR/VI SC* ABS ( EB ) ** ( 1 . -PLI ) * 

1  PRR+2.*EB*PRR 

YPRI ME ( 3 ) =-4 . *GSHEAR*EB-GSHEAR/VI SC*ABS ( EB ) ** ( 1 . -PLI ) * 
1  PZZ+-4.*EB*PZ2 

YPRIME  U ) =-2 . *EB*H 
YPRIME(5)=EB*R 

• 

RETURN 

END 

Q************* ************************** ******** 

SUBROUTINE  FCNJ(N,T, Y,PD) 

INTEGER  N 

REAL  Y(5) ,PD(N,N) ,T 

RETURN 

END 
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Q***************** **********************************  ****** 

C 

C  PROGRAM  PWLSLP.FOR 

C 

C  LUBRICATED  SQUEEZING  OF  POWER  LAW  FLUIDS. 

C 

Q****** ******************************************** ******* 

INTEGER  N,METH, MITER, INDEX, IWK(3),IER,IRV 

REAL  Y(3),WK(63),T,TOL,TEND,DELT 

DOUBLE  PRECISION  PRINT, THHO 

EXTERNAL  FCN,FCNJ 

COMMON/AB/A,B 

COMMON/AM/AMASS 

COMMON/G/GM 

COMMON/P/PLI 

Q******* ****************************** 

TYPE  10 

10  FORMATC  initial  CONDITIONS’/'  ENTER  TO,EBO ,H0 ,R0’ ) 
ACCEPT  *,T0,EB0,H0,R0 

TYPE  11 

11  FORMAT ('  ENTER  DENS,  VI SC  OR  CONSISTENCY  FACTOR,’ 

1  ’  POWER  LAW  INDEX,  MASS  OF  LOAD’) 

ACCEPT  * , DENS , VI SC , PLI , AMASS 
TYPE  12 

12  FORMAT (’  ENTER  TOLERANCE,  BELT’) 

ACCEPT  *,TOL,DELT 

TYPE  15 

15  FORMAT (’  ENTER  NPRT, INTERVAL’ ) 

ACCEPT  *,NPRT,AINTR 

TYPE  16 

16  FORMAT(’  ENTER  FILENAMES  FOR  PRINT  &  T  VS.  HHO ’ ) 
ACCEPT  11 7, PRINT, THHO 

117  FORMAT (AlO/AlO) 

^********* ************* ************** 

OPEN ( UNI T=  3 1 , DEVI CE= ’ DSK ’ , F I LE=PRI NT ) 

OPEN ( UNI T*  3  2 , DEVI CE* ’ DSK ’ , F I LE=THH0 ) 

IF{PLI .EQ.l. )  WRITE(31,17) 

IF(PLI .NE.l. )  WRITE(31,18) 

17  FORMAT(/’  LUBRICATED  SQUEEZING  OF  NEWTONIAN  FLUID’//) 

18  FORMAT(/’  LUBRICATED  SQUEEZING  OF  POWER  LAW  FLUID’//) 
IF(PLI .EQ.l. )  WRITE(31,25)  DENS , VI SC , AMASS 

IF(PLI .NE.l. )  WRITE(31,26)  DENS, VI SC, PL I , AMASS 

25  FORMAT (’  DENSITY  =  ’,E14.5/’  VISCOSITY  *  ’,E14.5/ 

1  '  MASS  =  ’ ,E14.5) 

26  FORMAT(’  DENS  =’,E14.5/'  CONSISTENCY  FACTOR  =',E14.5/ 

1  ’  POWER  LAW  INDEX  =  ’,E14.5/’  MASS  =  ’,E14.5) 

WRITE (31, 27)  TOL 

27  FORMATC  tolerance  =  ’,E14.5) 
0************************************ 

PI=3. 14159265357989 
N=3 
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METH=1 

MITER=0 

INDEX=1 

WRITE{31,28)  METH, MI TER, INDEX 
28  FORMAT{/'  METH  =  ',14/’  MITER  =  ',14/'  INDEX  =  ’,14) 


T=TO 

y(l)=EBO 
'  y(2)=H0 
Y(3)=R0 

GM=980.*AMASS 

A=PI*DENS 

B=PI*VISC 

EB=Y(1) 

H=Y(2) 

R-Y(3) 

HHO=H/HO 
WRITE (5, 20) 

WRITE(31,20) 

20  FORMAT (//7X, 'TIME' ,12X, 'EB' ,12X, 'H' ,13X, 'R' ,10X, 

1  'H/HO' ,10X, 'DELT'/) 

WRITE(5,21)  T,EB,H,R,HHO,DELT 
WRITE{31,21)  T,EB,H,R,HHO,DELT 
WRITE{32,*)  T,HH0 

21  FORMAT {6E14. 5) 

Q****************:iif*  *********  ********* 

DO  100  I=1,NPRT 
TEND=FLOAT(l )*AINTR 

CALL  DGEAR ( N , FCN , FCN J , T , DELT , Y , TEND , TOL , METH , MI TER , 
1  INDEX, IWK,WK,IER) 

EB*Y(1) 

H=Y(2) 

R=Y(3) 

HH  0 — H  /H  0 

WRI TE ( 3 1 , 2 1 )  TEND , EB , H , R , HH 0 , DELT 
WRITE(32,*)  TEND,HH0 
100  WRITE(5,21)  TEND, EB,H,R,HH0, DELT 

CLOSE ( UNI T=31) 

CLOSE ( UNI T*32) 

STOP 

END 

Q** ************************************ ******** 

SUBROUTINE  FCN(N,T, Y, YPRIME) 

C* ************************************** ''***** 

INTEGER  N 

REAL  Y(3) ,YPRIME(3) ,T 
COMMON /Ab/a,B 
COMMON/AM/AMASS 
COMMON/G/GM 
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COMMON/P/PLI 

Q******* ********* ********************** 

EB=y(i) 

EB2=EB*EB 

H=Y(2) 

H2=H*H 

R=Y(3) 

R2=R*R 
R4=R2*R2 
AAA=1 . 

IF(ABS(EB) .EQ.O. )  AAA=0. 

YPRIME ( 1 ) = ( GM+4 . *AMASS*H*EB2- A*R4*EB2/4 . 

1  +4.*A*R2*H2*EB2/3.-6.*B*R2* 

1  (12.)**((PLI-1.)/2.)*EB* 

1  ABS(EB)**(AAA*{PLI-1. ) ) ) 

1  /( A*R4/4 . +2 . *A*R2*H2/3 . +2 . *AMASS*H ) 

YPRIME(2)=-2.*EB*H 
YPRIME(3)«EB*R 

RETURN 

END 

(;*********************************************** 
SUBROUTINE  FCNJ(N,T, Y,PD) 

INTEGER  N 

REAL  Y(3)  ,PD(N,N)’,T 

RETURN 

END 
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Q******************************************************* 

C 

C  PROGRAM  JSSLP.FOR 
C 

c  LUBRICATED  COMPRESSIVE  FLOW  OF 

C  JOHNSON- SEGALMAN  MODEL. 

C 

Q******* ************** ************************ ********** 

INTEGER  N, METH, MITER, INDEX, IWK(5) ,IER,IRV 

REAL  Y(5) ,WK(115) ,T,TOL, TEND, DELT, MASS 

DOUBLE  PRECISION  PRINT, THHO 

EXTERNAL  FCN,FCNJ 

COMMON/AB/A,B 

COMMON/M/MASS 

COMMON/G/GM 

COMMON/VR/VI SC , RTIMEl , RTIME2 
COMMON/MP/G , RTIME , XI , EPS 

C************************************* 

7  TYPE  10 

10  FORMATC  INITIAL  CONDITIONS’/'  ENTER  T0,EB0,PRR0,PZZ0,H0,R0’ ) 
ACCEPT  *,TO,EBO,PRRO,PZZO,HO,RO 

TYPE  11 

11  FORMAT (’  ENTER  DENSITY,  AND  MASS  OF  LOAD') 

ACCEPT  *, DENS, MASS 

TYPE  111 

111  FORMATC  ENTER  SHEAR  MODULUS,  RELAXATION  TIME,  AND', 

1  '  RETARDATION  TIME' ) 

ACCEPT  *,G, RTIMEl, RTIME2 
TYPE  112 

112  FORMATC  ENTER  XI') 

ACCEPT  *,XI 

TYPE  12 

1 2  FORMAT ( ’  ENTER  TOLERANCE ,  DELT ' ) 

ACCEPT  *,TOL,DELT 

TYPE  15 

15  FORMATC  ENTER  NPRT,  INTERVAL' ) 

ACCEPT  *,NPRT,AINTR 

TYPE  16 

16  FORMATC  ENTER  FILENAMES  FOR  PRINT  /  T  VS.  HHO’) 

ACCEPT  117, PRINT, THHO 

117  FORMAT  (AIOAIO) 

*********************** ************ 

OPEN(UNIT=31,DEVICE=’DSK' ,FILB=PRINT) 
OPEN(UNIT=32,DEVICE='DSK' ,FILE=THH0) 

WRI TE (31  17) 

17  FORMAT(/’  LUBRICATED  COMPRESSIVE  FLOW  OF  JOHNSON-', 

1  'SEGALMAN  MODEL'//) 

WRITE{31,25)  DENS,G,RTIME1,RTIME2,XI ,MASS 
25  FORMATC  DENSITY  =  ',E14.5/'  SHEAR  MODULUS  =  ’, 

1  E14.5/'  RELAXATION  TIME  =  ',E14.5/ 

1  '  RETARDATION  TIME  =  ',E14.5/'  XI  =  ’ , 
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1  E14.5/'  MASS  =  ’ ,E14.5) 

WRITE (31, 27)  TOL 

27  FORMATC  tolerance  =  ',E14.5) 

Q********* ******************* *********** 

PI=3. 14159265357989 
N=5 

METH=2 

MITER=0 

INDEX=1 

WRITE(31,28)  METH, MI TER, INDEX 

28  FORMAT(/’  MWTH  =  ’,14/’  MITER  =  ’,14/’  INDEX  =  ’,14/) 
T=T0 

Y(1)=EB0 

Y(2)=PRR0 

Y{3)=PZZ0 

Y(4)=H0 

Y(5)=R0 

GM=980.*MASS 

A=PI*DENS 

EB=Y(1) 

PRR=Y(2) 

PZZ=Y(3) 

H*Y(4) 

R=Y(5) 

HH0=H/H0 

WRITE(5,20) 

WRITE(31,20) 

20  FORMAT (//7X, ’TIME’ ,12X, ’EB’ ,11X, ' PRR' ,11X, ' PZZ* , 

1  13X,’H' ,13X, ’R’ ,10X, ’H/HO’ ,10X, ’DELT’/) 

WRITE(5,21)  T, EB, PRR, PZZ, H,R,HH0, DELT 
WRITE(31,21)  T, EB, PRR, PZZ, H,R,HH0, DELT 
WRITE(32,*)  T,HH0 

21  FORMAT(8E14.5) 

Q************************************ 

DO  100  I=1,NPRT 
TEND=FLOAT( I ) *AINTR+T0 

CALL  DGEAR ( N , FCN , FCN J , T , DELT , Y , TEND , TOL , METH , MI TER , 

1  INDEX, IWK,WK,IER) 

EB=Y(1) 

PRR=Y(2) 

PZZ=Y(3) 

H=Y(4) 

R=Y(5) 

HHO=H/HO 

WRITE(31,21)  TEND, EB, PRR, PZZ, H,R,HH0, DELT 
WRITE (32,*)  TEND,HH0 

100  WRITE(5,21)  TEND, EB, PRR, PZZ, H,R,HH0, DELT 

CLOSE(UNIT=31,DISPOSE=’ PRINT’ ) 
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CLOSE(UNIT=32) 

STOP 

END 

Q********************************************** 

SUBROUTINE  PCN(N,T, Y, YPRIME) 

Q**** **************************** ************* 

INTEGER  N 

REAL  Y(5) ,YPRIME(5) ,T,MASS 
COMMON/AB/A,B 
COMMON/M/MASS 
COMMON/G/GM 

COMMON/VR/VI SC , RTIMEl , RTIME2 
COMMON/MP/G , RTIME , XI 

Q*********T*r****  *******:*:****  ************ 

PI=3. 14159265357989 
EB=Y(1) 

EB2=EB*EB 

PRR=Y(2) 

PZZ=Y(3) 

H=Y(4) 

H2=H*H 

R=Y(5) 

R2=R*R 

R4»R2*R2 

YPRIME ( 1 ) = (GM+4 . *MASS*H*EB2-A*R4*EB2/4 . 

1  +4.*A*R2*H2*EB2/3.-PI*R2*(PRR-PZZ) ) 

1  /( A*R4/4 . +2 . *A*R2*H2/3 . +2 . *MASS*H) 

YPRIME { 2 ) =2 . *G* ( EB+RTIME2* ( YPRIME ( 1 ) -2 . * ( 1 . -XI ) 

1  *EB2 ) ) +2 . * ( 1 . -XI ) *EB*PRR-PRR/RTIME1 

YPRIME ( 3 ) =-4 . *G* ( EB+RTIME2* (YPRIME ( 1 ) +4 . * ( 1 . -XI ) 
1  *EB2 ) ) -4 . * ( 1 . -XI ) *EB*PZZ-PZZ/RTIME1 

YPRIME (4 )=-2.*EB*H 
YPRIME ( 5) =EB*R 

RETURN 

END 

Q*********** ******* ***************************** 

SUBROUTINE  FCNJ{N,T, Y,PD) 

Q** ***************************** **************** 

INTEGER  N 

REAL  Y(5) ,PD(N,N) ,T 

RETURN 

END 
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Q*** ************************************************ ********** 

C 

C  PROGRAM  NJSSLP.FOR 
C 

C  LUBRICATED  SQUEEZING  OF  INELASTIC  FLUIDS 

C  WITH  JOHNSON-SEGALMAN  VISCOSITY. 

C 

0************** ************** ********************************* 
INTEGER  N,METH,MITER, INDEX, IWK{3 ) ,IER,IRV 
REAL  y{3),WK(63),T,TOL,TEND,DELT 
DOUBLE  PRECISION  PRINT, THHO 
EXTERNAL  FCN,FCNJ 
COMMON/AB/A,B 
COMMON/AM/AMASS 
COMMON/G/GM 

COMMON/RX/RTIMEl , RTIME2 , XI 
0************************^** ********** 


10 

TYPE  10 

FORMATC  INITIAL  CONDITIONS'/'  ENTER 

T0,BB0,H0,R0' ) 

11 

ACCEPT  *,TO,EBO,HO,RO 

TYPE  11 

FORMATC  ENTER  DENS,  ZERO  VI SC,  MASS 

OF  LOAD' ) 

111 

ACCEPT  *, DENS, VI SCO, AMASS 

TYPE  111 

FORMATC  ENTER  RTIME1,RTIME2, XI '  ) 

12 

ACCEPT  *,RTIME1,RTIME2,XI 

TYPE  12 

FORMATC  ENTER  TOLERANCE,  DELT' ) 

15 

ACCEPT  *, TOL, DELT 

TYPE  15 

FORMATC  ENTER  NPRT,  INTERVAL' ) 

16 

ACCEPT  *,NPRT,AINTR 

TYPE  16 

FORMATC  ENTER  FILENAMES  FOR  PRINT  & 

T  VS .  HHO ' )  - 

117 

ACCEPT  11 7, PRINT, THHO 

FORMAT (AlO/AlO) 

0****************** ******* *********** 

OPEN ( UNI T= 31, DEVI CE='DSK’ ,FILE=PRINT) 

OPEN ( UNI T=  3  2 , DEVI CE= ’ DSK ' , FI LE=THH0 ) 

WRITE(31,17) 

17  FORMAT(/’  LUBRICATED  SQUEEZING  OF  INELASTIC  FLUIDS', 

1  /’  WITH  JOHNSON-SEGALMAN  VISCOSITY') 

WRITE(31,18)  DENS, VI SCO, AMASS 

18  FORMATC  DENSITY  =  ',E14.5/'  ZERO  VISCOSITY  =  ',E14.5/ 

1  '  MASS  OF  LOAD  =  ',E14.5) 

WRITE(31,25)  RTIME1,RTIME2,XI 
25  FORMATC  RELAXATION  TIME  *  ',E14.5/ 

1  '  RETARDATION  TIME  =  ',E14.5/’  XI  =  ',E14.5) 

WRITE(31,27)  TOL 

27  FORMATC  TOLERANCE  =  ',E14.5) 

0************* ********* ************** 


169 


PI=3. 14159265357989 
N=3 

METH=1 

MITER=0 

INDEX=1 

,,  -  ■  -»/• 

T=TO 

Y(1)=EB0 

y(2)=H0 

y(3)=R0 

GM=980.* AMASS 

A=PI*DENS 

B=PI*VISCO 

EB=y(l) 

H=y(2) 

R=Y(3) 

HHO=H/HO 
WRITE (5, 20) 

WRITE(31,20)  fool  1 -5V  »H'  1  'R'  lOX, 

20  F0RMAT(//7X,’TIMEM2X,  EB  ,12X,  H  ,13X,  R  ,1U  , 

1  ’H/HO’ ,10X,’DELT’/) 

WRITE{5,21)  T,EB,H,R,HH0,DELT 
WRITE(31,21)  T,EB,H,R,HH0 ,DELT 
WRITE(32,*)  T,HH0 

21  F0RMAT(6E14.5) 

C*********************************  * 

DO  100  I=l,NPRT 

C ALl'dGEAR ( N ! FCN^FCN J , T , DELT , Y , TEND , TOL , METH , MI TER , 
1  INDEX, IWK,WK,IER) 

EB=Y(1) 

H=Y(2) 

R=Y(3) 

hho=h/ho 

WRITE(31,21)  TEND, EB,H,R,HH0, DELT 
WRITE(32,*)  TEND,HH0 
100  WRITE(5,21)  tend, EB,H,R,HH0, DELT 

CLOSE(UNIT=31,DISPOSE=' PRINT  ) 

CLOSE (UNIT=32) 

STOP 

c.....5«*>**********«******‘****************** 

SUBROUTINE  FCN{N,T,Y,YPRIME) 

c************************************ 

INTEGER  N  ' 

REAL  Y(3),YPRIME(3),T 
COMMON/AB/A,B 
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COMMON/AM/AMASS 
COMMON/G/GM 

COMMON/RX/RTIMEl , RTIME2 , XI 

Q****** *************** ***************** 

EB=y(i) 

EB2»EB*EB 
H=Y(2) 

H2=H*H 

R=y(3) 

R2=R*R 
R4=R2*R2 

yPRIME(l) 

1 
1 
1 
1 
1 
1 

yPRIME(2) 
yPRIMEO) 

RETURN 

END 

Q********* ************************** ************ 

SUBROUTINE  FCNJ(N,T,y ,PD) 

Q*** ************************************** ****** 

INTEGER  N 

REAL  y(3) ,PD(N,N) ,T 

RETURN 

END 


= {GM+4 . *AMASS*H*EB2-A*R4*EB2/4 . 

+4 . *A*R2*H2*EB2/3 . -6 . *B*R2*EB 
* ( 1 . +12 . *XI * ( 2 . -XI ) *RTIME1* 
RTIME2*EB2) 

/(l.+12.*XI*(2.-XI )*RTIME1* 
RTIME1*EB2)) 

/( A*R4/4 . +2 , *A*R2*H2/3 . +2 . *AMASS*H) 
=-2.*EB*H 
=EB*R 
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c*********************************************************** 

C 

C  PROGRAM  SMMSLP.FOR 

C 

c  LUBRICATED  COMPRESSIVE  FLOW  OF 

C  MARRUCCI  STRUCTURAL  MODEL. 

C 

c******* *********************************** ***************** 
INTEGER  N,METH,MITER,INDEX,IWK{6) ,IER,IRV 
REAL  Y{6) ,WK(144) ,T,TOL, TEND, DELT, MASS 
DOUBLE  PRECISION  PRINT, THHO 
EXTERNAL  FCN,FCNJ 
COMMON/AB/A,B 
COMMON/M/MASS 
COMMON/G/GM 

COMMON/VR/VI SC , RTIMEl , RTIME2 
COMMON/MP/GO , RTIMEO , XI , EPS 
(;************************************* 

TYPE  1 

I  FORMATC  LUBRICATED  COMPRESSIVE  FLOW  OF'/ 

1  ’  MARRUCCI  STRUCTURAL  MODEL'// 

1  '  ALL  INPUT  DATA  IN  CGS  UNITS'/) 

7  TYPE  10 

10  FORMATC  INITIAL  CONDITIONS'/ 

1  '  ENTER  TO, EB0,PRR0,PZZ0, H0,R0,X0' ) 

ACCEPT  * , TO , EBO , PRRO ,PZZ0,H0,R0,X0 
TYPE  11 

II  FORMATC  ENTER  DENSITY,  AND  MASS  OF  LOAD') 

ACCEPT  *, DENS, MASS 

TYPE  111 

III  FORMATC  ENTER  ZERO  SHEAR  MODULUS,  ZERO  RELAX  TIME') 
ACCEPT  *, GO, RTIMEO 

TYPE  12 

12  FORMATC  ENTER  TOLERANCE,  DELT') 

ACCEPT  *,TOL,DELT 
TYPE  15 

1 5  FORMAT { '  ENTER  NPRT , I NTERVAL ' ) 

ACCEPT  *,NPRT,AINTR 

TYPE  16 

16  FORMATC  ENTER  FILENAMES  FOR  PRINT  /  T  VS.  HHO ' ) 

ACCEPT  11 7, PRINT, THHO 

117  FORMAT (AIO/AIO) 

(;*****  *************mr  ***********  ****** 

OPEN ( UNI T=  3 1 , DEVI CE= ' DSK ' , F I LE=PRI NT ) 
OPEN(UNIT=32,DEVICE='DSK' ,FILE=THH0) 

WRITE(31,17) 

17  FORMATC/'  LUBRICATED  COMPRESSIVE  FLOW  OF  ', 

1  'MARRUCCI  STRUCTURAL  MODEL'//) 

WRITE(31,25)  DENS, GO, RTIMEO, MASS 
25  FORMATC  DENSITY  =  ',E14.5/'  ZERO  SHEAR  MODULUS  =  ', 

1  E14.5/'  ZERO  RELAXATION  TIME  =  ',E14.5/ 
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1  '  MASS  =  ' ,E14.5) 

WRITE (31, 27)  TOL 

27  FORMATC  tolerance  =  •,E14.5) 

Q* ************************************** 

PI=3. 14159265357989 
N*6 

METH=2 

MITER=0 

INDEX*1 

WRITE(31,28)  METH, MI TER, INDEX 

28  FORMATC/'  METH  =  ',14/'  MITER  =  ',14/'  INDEX  =  ’,14/) 
T=T0 

Y(1)=EB0 

Y(2)=PRR0/(G0*X0) 

Y(3)=PZZO/(GO*XO) 

Y(4)=X0 

Y(5)=H0 

Y{6)=R0 

GM=980.*MASS 

A=PI*DENS 

EB=Y(1) 

PRRH=Y(2) 

PZZH»Y(3) 

X=Y(4) 

H=Y(5) 

R=Y(6) 

G=G0*X 
PRR=PRRH*G 
PZZ*PZZH*G 
HH0=H/H0 
WRITE (5, 20) 

WRITE (31, 20) 

20  FORMAT (//7X, 'TIME' ,12X,'EB’ ,11X,’PRR' ,11X,'PZZ' , 

1  13X, 'H' ,13X, 'R' ,10X, 'H/HO' ,10X, 'DELT' ,12X, 'X'/) 

WRITE(5,21)  T,EB,PRR,PZZ,H,R,HHO,DELT,X 
WRITE (31, 21)  T,EB,PRR,PZZ,H,R,HHO,DELT,X 
WRITE(32,*)  T,HH0 

21  FORMAT(9E14.5) 

C************************************ 

DO  100  I=1,NPRT 
TEND=FLOATh  )*AINTR+T0 

CALL  DGEAR ( N , FCN , FCN J , T , DELT , Y , TEND , TOL , METH , MI TER , 

1  INDEX, IWK,WK,IER) 

EB=Y(1) 

PRRH=Y(2) 

PZZH*Y(3) 

X=Y(4) 

H=Y(5) 
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R=Y(6) 

G=GO*X 

PRR*PRRH*G 

PZZ=PZZH*G 

hho=h/ho 

WRITE (31, 21)  TEND, EB,PRR,PZZ,H,R,HHO, BELT, X 
WRITE (32,*)  TEND,HHO 

100  WRITE(5,21)  TEND, EB,PRR,PZZ,H,R,HHO, BELT, X 

CLOSE (UNIT=31) 

CLOSE (DNIT=32) 

STOP 

END 

Q***ii^*  *********************************  ******** 

SUBROUTINE  FCN(N,T,Y,y PRIME) 

Q*************** *********************** ******* 

I NTEGER  N 

REAL  Y(6) ,YPRIME(6) ,T,MASS 
COMMON/AB/A,B 
COMMON/M/MASS 
COMMON/G/GM 

COMMON/VR/VI SC , RTIMEl , RTIME2 
COMMON/MP/GO , RTIMEO , XI , EPS 

Q* ************************************* 

PI=3. 14159265357989 
EB=y(l) 

EB2»EB*EB 

PRRH=Y(2) 

PZZH=y(3) 

X=Y(4) 

X14=X**1 • 4 

RX14=RTIME0*X14 

H=y(5) 

H2»H*H 

R=Y(6) 

R2=R*R 

R4=R2*R2 

yPRIME ( 1 ) = (GM+4 . *MASS*H*EB2-A*R4*EB2/4 . 

1  +4.*A*R2*H2*EB2/3.-PI*R2*G0*X*(PRRH-PZZH) ) 

1  /( A*R4/4 . +2 . *A*R2*H2/3 . +2 . *MASS*H) 

YPRIME ( 2 ) =2 . *EB-PRRH/RX14+2 . *EB*PRRH 

YPRIME( 3 ) =-4 . *EB-PZZH/RX14-4 . *EB*PZZH 

YPRIME ( 4 ) = ( 1 . -X) /RX14-0 . 25*X*SQRT ( 2 . *PRRH+PZZH) /RX14 

yPRIME(5)=-2.*EB*H 

yPRIME(6)=EB*R 

RETURN 

END 

(;*********************************************** 

SUBROUTINE  FCNJ(N,T,Y,PD) 

INTEGER  N 
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REAL  Y(6) ,PD(N,N) ,T 

RETURN 

END 
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Q*********************************************************** 

c 

C  PROGRAM  NSMSLP.FOR 
C 

C  LUBRICATED  SQUEEZING  OF  INELASTIC  FLUIDS 

C  WITH  STRUCTURE-DEPENDENT  VISCOSITY. 

C 

C************** ********************************************* 

INTEGER  N,METH,MITER, INDEX, IWK(3) ,IER,IRV 

REAL  Y(3),WK(63),T,TOL,TEND,DELT 

DOUBLE  PRECISION  PRINT, THHO 

EXTERNAL  FCN,FCNJ 

COMMON/AB/A,B 

COMMON/AM/AMASS 

COMMON/G/GM 

COMMON/RT/RTI MEO 

COMMON/X/X 

C****** ********************** ********* 

TYPE  10 

10  FORMATC  initial  CONDITIONS'/*  ENTER  TO ,EB0 ,H0 ,R0 ,X0 ' > 
ACCEPT  *,T0,EB0,H0,R0,X0 

TYPE  11 

11  FORMATC  ENTER  DENSITY,  MASS  OF  LOAD') 

ACCEPT  *, DENS, AMASS 

TYPE  111 

111  FORMATC  ENTER  ZERO  VISCOSITY,  ZERO  RELAXATION  TIME’) 
ACCEPT  *,  VI SCO,  RTIMEO 
TYPE  12 

12  FORMATC  ENTER  TOLERANCE,  DELT' ) 

ACCEPT  *,TOL,DELT 

TYPE  15 

15  FORMATC  ENTER  NPRT,  INTERVAL' ) 

ACCEPT  *,NPRT,AINTR 

TYPE  16 

16  FORMATC  ENTER  FILENAMES  FOR  PRINT  &  T  VS.  HHO') 

ACCEPT  11 7, PRINT, THHO 

117  FORMAT (AlO/AlO) 

Q* ****** ***************************** 

OPEN ( UNI T= 31, DEVI CE»’DSK' ,FILE*PRINT) 

OPEN ( UNI T= 3  2 , DEVI CE= ' DSK ’ , FI LE=THH0 ) 

WRITE(31,25) 

25  FORMAT(/'  LUBRICATED  SQUEEZING  OF  INELASTIC  FLUIDS' 

1  /’WITH  STRUCTURE -DEPENDENT  VISCOSITY.') 

WRI TE ( 31 , 26 )  DENS , VI SCO , AMASS , RTIMEO 

26  FORMAT(/'  DENS  =  ',E14.5/’  ZERO  VI SC  =  ',E14.5/ 

1  ’  MASS  OF  LOAD  =  ',E14.5/ 

1  '  ZERO  RELAXATION  TIME  =  ',E14.5) 

WRITE (31, 27)  TOL 

27  FORMATC  TOLERANCE  =  ',E14.5) 

Q********************* *************** 

PI=3. 14159265357989 
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N=3 

METH=1 

MITER=0 

INDEX=1 

WRITE(31,28)  METH, MI TER, INDEX 
28  FORMAT(/’  METH  =  ’,14/'  MITER  =  ',14/'  INDEX  =  ',14) 

T=TO 

y(l)=EBO 

Y(2)=H0 

Y(3)»R0 

X=XO 

GM=980.*AMASS 

A=PI*DENS 

B=PI*VISCO 

EB=Y(1) 

H=Y(2) 

R=Y(3) 

HHO=H/HO 
WRITE (5, 20) 

WRITE (31, 20) 

20  FORMAT {//7X, 'TIME' ,12X, 'EB' ,12X, 'H' ,13X, 'R' ,10X, 

1  'H/HO' ,10X, 'DELT' ,12X, 'X'/) 

WRITE(5,21)  T,EB,H,R,HHO,DELT,X 
WRITE(31,21)  T,EB,H,R,HHO,DELT,X 
WRITE(32,*)  T,HH0 

21  FORMAT (7E14. 5) 

Q* ***************************  ******** 

DO  100  I=1,NPRT 
TEND»FLOAT( I ) *AINTR 

CALL  DGEAR ( N , FCN , FCN J , T , DELT , Y , TEND , TOL , METH , MI TER , 

1  INDEX, IWK,WK,IER) 

EB=Y{1) 

H=Y(2) 

R=Y(3) 

HHO=H/H0 

WRITE (31, 21)  TEND, EB,H,R,HH0, DELT, X 
WRITE (32,*)  TEND,HH0 

100  WRITE(5,21)  TEND, EB,H,R,HH0, DELT, X 

CLOSE ( UNI T= 31) 

CLOSE (UNI T=3 2) 

STOP 

END 

^******************************* *************** 

SUBROUTINE  FCN(N,T, Y,YPRIME) 

^********************** ****** ***************** 

INTEGER  N 

REAL  Y(3) ,YPRIME(3) ,T 
COMMON/AB/A , B 
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COMMON/AM/ AMAS  S 
COMMON/G/GM 
COMMON /RT/RTI MEO 

c************************************** 

EB=Y(1) 

EB2=EB*EB 

H=Y(2) 

H2=H*H 

R=Y(3) 

R2=R*R 

1002  ^=X-(i!38564*RTIME0*EB*X**2.4+X-1.)/ 

lUU^  iN  A  ,2.4*1.38564*RTIME0*EB*X**1.4+1.) 

IF(ABS( (XN-X)/XN) .LT. 0.001)  GO  TO  1001 

1005  X=XN 

GO  TO  1002 
1001  X=XN 

1000  YPRIME  ( 1 )  =  {GM+4 .  S*H*EB2-A*R4*EB2/4 . 

1  +4.*A*R2*H2*EB2/3.-6.*B*R2* 

1  (1  -X)/(1.38564*RTIME0) )  ^  ^ 

1  /( A*R4/4 . +2 . *A*R2*H2/3 . +2 . *AMASS*H) 

YPRIME {2) --2. *EB*H 
YPRIME(3)=EB*R 

RETURN 

c.....?S.......***»***********‘**************** 

INTEGER  N 

REAL  Y(3),PD(N,N),T 

RETURN 

END 


APPSI^DIX  P 


RHEOLOGICAL  DATA 


F.l  Oscillatory  shear  data  on  silicone  polymer 


(1)  Temperature  =  18°C 


CIRCULAR 

STORAGE 

LOSS 

COMPLEX 

FREQUENCY, <0 

MODULUS, g' 

MODULUS,  g'''’ 

VISCOSITY, 7* 

(rad/ SEC) 

(DYN/CM2) 

(DYN/CM2) 

(POISE) 

0.1 

0.750E+04 

0.113E+07 

0.113E+07 

0.31 

0.400E+05 

0.360E+07 

0.117E+07 

1.0 

0.275E+06 

0.107E+07 

O.llOE+07 

3.1 

0.145E+07 

0.230E+07 

0.877E+06 

10. 

0.395E+07 

0.275E+06 

0.481E+06 

31. 

0.600E+07 

0.185E+06 

0.203E+06 

100. 

0.620S+07 

0.860E+06 

0.626E+05 

Horizontal 

shift  factor  (^j) 

=  1.60 

Temperature 

=  23°C 

CIRCULAR 

STORAGE 

LOSS 

COMPLEX 

FREQUENCY, 

MODULUS,  g' 

MODULUS,  g'" 

VISCOSITY, 7 

(RAD/ SEC) 

(DYN/CM2) 

(DYN/CM2) 

(POISE) 

0.1 

0.240E+04 

0. 730E+05 

0.730E+06 

0.31 

0.130E+05 

0.230E+06 

0.743E+06 

1.0 

O.lOlE+06 

0.700E+06 

0.707E+06 

3.1 

0.710E+06 

0.185E+07 

0.639E+06 

10. 

0.296B+07 

0.310E+07 

0.429E+0O 

31. 

0.590E+07 

0.265E+07 

0.208E+06 

100. 

0.650E+07 

0.135E+07 

0.564E+05 

Horizontal 

shift  factor  (<^t) 

=  1.0 
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Temperature 

=  3  5°C 

CIRCULAR 
FREQUENCY,^ 
(RAD/ SEC) 

STORAGE  , 

MODULUS , G 
(DYN/CM2) 

LOSS  ,, 

MODULUS , G 
(DYN/CM2) 

COMPLEX  ^ 

VISCOSITY, 

(POISE) 

0.1 

0.31 

1.0 

3.1 

10. 

31. 

100. 

O.lOOE+03 

0.305E+04 

0.109E+05 

0.173S+05 

0.115S+07 

0. 390E+07 
0.570E+07 

0.325E+05 

O.lOOE+06 

0.330S+06 

0.940E+06 

0.230E+07 

0.320E+07 

0.210E+07 

0.325E+06 

0.323E+06 

0.330E+06 

0.308E+06 

0.257E+06 

0.163E+06 

0.607S+05 

Horizontal 

shift  factor 

=  0.44 
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F. 2  Steady  shear  data  on  TLA-227 


Temperature 

=  23. 5*^0 

- 

SHEAR 

SHEAR 

NORMAL 

REMARK 

RATE,  t 

STRESS,  7.a 

STRESS,  7»-7tx 

(SEC') 

(DY£J/CM^  ) 

(DYN/CM^) 

32.1 

0.934E+04 

CAPILLARY 

37.6 

0. 105E+05 

- 

(  L/r  s  s.68  ) 

51.0 

0.136E+05 

- 

it 

74.8 

0. 192E+05 

- 

// 

118.2 

0.286E+05 

- 

if 

Temperature 

SHEAR 

=  2  7®C 

SHEAR 

NORMAL 

REMARK 

RATE,  t 

STRESS, 

STRESS, 

(SEC“') 

(dyn/cm"^  ) 

(DYN/CM^) 

1.0 

0.260E+03 

O.llOE+03 

CONE  &  PLATE 

2.5 

0.600E+03 

0. 300E+03 

// 

6.3 

0.145E+04 

0.120E+04 

It 

16. 

0.336E+04 

0.440E+04 

U 

25. 

0.463E+04 

0.880E+04 

// 

34.5 

0. 751E+04 

- 

CAPILLARY 

70.5 

0.140S+05 

- 

93. 

1.83E+05 

- 

// 

122.8 

6.224E+05 

- 

it 

168. 

0. 291E+05 

- 

it 

195.3 

0.333E+05 

- 

It 

62.8 

0.119E+05 

- 

CAPILLARY 

123.3 

0.212E+05 

- 

(L/r  =  /60) 

194.5 

0. 316E+05 

- 

// 

259.1 

0.423E+05 

- 

It 

357.9 

0. 557E+05 

- 

// 
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Temperature 

=  30®C 

SHEiAR 

RATE,  7 
(3EG"') 

SHEAR 

STRESS ,  Tiz. 
(DYN/CM^) 

NORMAL 

STRESS,  7ii-7z2. 
(DYN/GM'^) 

REMARK 

1.6 

2.5 

6.3 

10. 

16 . 

25. 

0.304E+03 

0.500S+03 

0.113E+04 

0.185S+04 

0.272E+04 

0.400E+04 

0.800E+02 

0.180E+03 

0.150E+04 

0.320S+04 

0.533S+04 

CONE  Sc 

ft 

t* 

ft 

it 

ft 

PLATE 

)  Temperature 

=  35®G 

SHEAR  . 

RATE, 

(SEC'*) 

SHEAR' 

STRESS,  7iz 
(DYH/GM^) 

NORMAL 
STRESS,'?',, -7m. 

(dyn/gm-^) 

REMARK 

1.6 

4. 

10. 

0.216E+03 

0.4838+03 

0.115E+04 

0.206S+03 

0.940S+03 

CONE  Sc 

// 

ft 

PLATS 
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F.3  Steady  shear  data  on  3 ♦ 3 


Temperature 

=  25  °  C 

SHEAR 

SHEAR 

RATE,t 

STRESS,  Ta 

( 3EC’' ) 

(DYN/CM*) 

1.05 

0.262E+03 

2.1 

0. 330E+03 

5.28 

0.447E+03 

10.5 

0.561E+03 

21. 

0.683E+03 

52.8 

0. 879E+03 

105. 

0.105S+04 

210. 

0. 126E+04 

419. 

0.151E+04 

wt.  %  PAA  in  water  solution 


NORMAL  REMARK 

STRSSS,7;.-?2i 

(DYN/CM^) 

CONE  &  PLATE 


—  u 

~  » 

~  tt 

0.574E+04 

0.362E+04 

0.137E+05  '' 

0.211S+05  » 


APPENDIX  G 


EXPERIMENTAL  DATA  ON  THE  SQUEEZING  FLOW 
G • 1  Newtonian  Fluids 

(1)  Viscasil  SOOOO 

Run  1  : 


R  = 

1.252  cm 

2Ho 

=  0.102  cm 

m  = 

7  = 

6239  g 

530  poise 

(at  25®C) 

t  (sec) 

2H  ( cm) 

0. 

0.1021 

0.0333 

0.0760 

0.102 

0.0534 

0.204 

0.0455 

0.305 

0.0335 

0.403 

0.0340 

0.612 

0.0232 

0.316 

0.0247 

Run  2  : 

R  =  1.252  cm 
2Ho  =  0.101  cm 
m  =  6239  g 

=  530  poise  (at  25°G) 


t  (sec) 

2H  ( cm) 

0. 

0.101 

0.0403 

0.0750 

0.112 

0.0569 

0.214 

0.0443 

0.316 

0.0375 

0.413 

0.0329 

0.520 

0.0297 
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(2)  Dow  Corning  200  fluid,  12500 


Run  1 


R  =  2 . 534  cm 
2Ho  =  0.211  cm 
m  =  3676  g 

^  =  130  poise  (at  22.5°C) 


t  ( sec) 

2H  (cm) 

0. 

0.211 

0.0214 

0.138 

0.0897 

0.148 

0.2254 

0.112 

0.3274 

0.0969 

0.429 

0.0855 

0.5314 

0.0775 

0.633 

0.0719 

Run  2 


R  =  1.252  cm 

2Ho  =  0.122  cm 

m  =  6290  g 

'tj  =  121  poise  (at 

25.5°G) 

t  (sec) 

2H  (cm) 

0. 

0.122 

0.0612 

0.0403 

0.1632 

0.0244 

0.2652 

0.0201 

0.3672 

0.0171 
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G.2  Viscoelastic  materials 


(1)  Silicone  polymer 


Run  1 


R  =  0. 505  cm 


2Ho  =  0.123  cm 
m  =  17480  g 
^  =  920000  poise 

(at  20®C) 

X  =  0.12  sec 

• 

t  (sec) 

2H  (cm) 

0. 

0.123 

0.02 

0.111 

0.05 

0.109 

6.25 

0.0986 

0.55 

0.0358 

0.95 

0.0740 

1.55 

0.0617 

2.05 

0.0543 

2.55 

0.0486 

3.05 

0.0444 

Vmftx  =  1.  cm/ sec 
TSma,K  ~  105.  sec^ 

t’maxA  ~  12.6 
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Run  2 


R  =  0 . 505  cm 
2Ho  =  0.26  cm 
m  =  8612  g 

^  =  655000  poise  (at  24.4°C) 
X  =  0 . 1  sec 


t  (msec) 

0. 

9.2 

14.3 

18.9 

21.9 

24.5 
27.8 

31.1 

34.4 

35.7 

39.8 

44.9 

47.4 
50.0 

52.5 

55.6 
61.0 

67.9 

75.5 

80.6 

84.7 

90.3 

93.8 

102. 

111.2 

116.3 

134.1 
152. 

162.2 

172.9 


2H  (cm) 

0.26 
0.258 
0.250 
0.240 
0.230 
0.224 
0.220 
0.224 
0.230 
0.234 
0.236 
0.234 
0.  230 
0.224 
0.220 
0.218 
0.220 
0.224 
0.220 
0.216 
0.215 
0.216 
0.216 
0.214 
0.211 
0.211 
0.207 
0.204 
0.203 
0.20 


Vmttx  =  3.3  cm/  sec 
2^in«uc  —  90.  sec  ^ 
j  m«.K  X  =  9  . 
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(2)  TLA-227 


Run  1 


R  =  2.534  cm 
2Ho  =  0.14  cm 
ra  =  9307  g 

^  =  455(t)*’^*"'  (at  23.4°C) 


t  (sec) 

2H  (cm) 

0. 

0.14 

0.0306 

0.126 

0.102 

0.1124 

0.161 

0.1005 

0.306 

0.0864 

0.510 

0.0724 

0.714 

0.0636 

0.913 

0.0570 

1.122 

0.0519 

V  maoc~  0 . 46 

cm/ sec 

218. 

sec”* 

^i*wwA“ 


Run  2 

R  =  1.25  cm 


2Ho  =  0.0584  cm 

m  =  3143  g  , 

^  =  455(7)®'®^”' 

t  (sec) 

2H  (cm) 

0. 

6.0584 

0.01275 

0.0549 

0.06375 

0.0490 

0.1148 

0.0458 

0.1658 

0.0429 

0.2168 

0.0404 

0.2678 

0.0382 

Vmoix  =  0.275  cm/ sec 
=  347  .  sec"! 

1  =2.95 
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Run  3 


R  =  1.252  cm 
2Ho  =  0.143  cm 
m  =  3143  g 
^  ^  455  it)"- 


t  (sec) 

2H  (cm) 

0. 

0.148 

0.02 

0.105 

0.025 

0. 103 

0.032 

0.0986 

0.065 

0.0829 

0.0875 

0.074 

V|ti*x  ~  2.16 

cm/ sec 

• 

7w<xx  ~  5  94 . 
irmiwA  =3.6 

sec"* 

189 


(3)  PAA-v/ater  solution 


Run  1 


R  =  6.35  cm 
2Ho  =  0.2915  cm 
m  =  4751  g 

=  302.7(jr)°-^^^  “* 


t  (sec) 

2H  (cm) 

0. 

0.2915 

0.0115 

0.2693 

0.0265 

0.261 

0.0335 

0.2584 

0.053 

0.2466 

0.112 

0.2145 

0.162 

0.1864 

0.212 

0.1618 

0.262 

0.1435 

0.312 

0.131 

Vrt,«x=  1-93 

cm/ sec 

Ifmox  ~  893 . 

8.7 

sec”* 

